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Abstract 

We analyze the symmetries and the self-consistent perturbative 
approaches of dynamical field theories for glass- forming liquids. In 
particular we focus on the time-reversal symmetry which is crucial 
to obtain fluctuation dissipation relations (FDR). Previous field theo- 
retical treatment violated this symmetry, whereas others pointed out 
that constructing symmetry preserving perturbation theories is a cru- 
cial and open issue. In this work we solve this problem and then 
apply our results to the Mode Coupling Theory of the glass transition 
(MCT). 

We show that in the context of dynamical field theories for glass- 
forming liquids time-reversal symmetry is expressed as a non-linear 
field transformation that leaves the action invariant. Because of this 
non-linearity standard perturbation theories generically do not pre- 
serve time-reversal symmetry and in particular fluctuation-dissipation 
relations. We show how one can cure this problem and set up sym- 
metry preserving perturbation theories by introducing some auxiliary 
fields. As an outcome we obtain Schwinger-Dyson dynamical equa- 
tions that preserve automatically FDR and that serve as a basis for 
carrying out symmetry preserving approximations. We apply our re- 
sults to the Mode-Coupling Theory of the glass transition revisiting 
previous field theory derivations of MCT equations and showing that 
they generically violate FDR. We obtain symmetry preserving mode- 
coupling equations and discuss their advantages and drawbacks. Fur- 
thermore we show, contrary to previous works, that the structure of 
the dynamic equations is such that the ideal glass transition is not 
cut-off at any finite order of perturbation theory, even in the presence 
of coupling between current and density. The opposite results found 
in previous field theoretical works, as the ones based on non-linear 
fluctuating hydrodynamics, were only due to an incorrect treatment 
of time-reversal symmetry. 
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Introduction 



Liquids, if cooled fast enough in order to avoid crystallization, enter gener- 
ically in a metastable phase in which remarkable dynamical phenomena oc- 
cur E] upon decreasing the temperature. For instance the structural re- 
laxation time increases very fast, in many cases faster than an Arrhenius law, 
correlation functions have slower than exponential relaxations, and dynamics 
become strongly heterogeneous, contrary to the simple motion which occurs 
at higher temperature. Typically at a temperature that is roughly two third 
of the melting temperature the relaxation time becomes macroscopic, of the 
order of minutes or hours, and the liquid freezes into an amorphous solid 
called glass. 

The equilibrium dynamics of moderately supercooled liquids is rather well 
described by the Mode-Coupling Theory |21 IH E] (MCT). MCT is a closure 
scheme for the correlation function developed twenty years ago [HI Ej lead- 
ing to a nonlinear integro-differential equation which has to be solved self- 
consistently. It predicts a power law divergence of the relaxation time and of 
the viscosity at a finite temperature Tmct- Although it is now clear that this 
is a spurious transition several quantitative results compare well with experi- 
mental and numerical findings as for example the wave-vector dependence of 
the Debye- Waller factor (the analog of the Edwards- Anderson parameter for 
spin-glasses) see IH E] • The conventional interpretation is that at T M ct 
there is a dynamical crossover: "hopping or activated" events not contained 
in the theory and which can be roughly neglected for T > Tmct prevent the 
existence of the transition and dominate the slowing down for T < Tmct- 
Thus, although the transition is strictly speaking avoided, the dynamics for 
T > Tmct can be explained in terms of MCT. This interpretation is based 
on theoretical extensions of MCT [HI |H] and recent developments originating 
from mean- field disordered systems [TO] - However, it is important to stress 
that the MCT dynamical crossover is not that sharp in real systems, therefore 
many MCT statements have to be considered with great care. Actually even 
the existence of a MCT crossover remains a long debated issue in glass litera- 
ture, see e.g. ^2]- MCT was derived originally using the projection operator 
formalism of Mori and Zwanzig [7] . A few years later, MCT was rederived 
starting from stochastic nonlinear hydrodynamics equations as a one-loop 
self-consistent theory jHj. This field theory derivation and subsequent ones 
have been criticized for two reasons. First because the mechanism behind the 
MCT transition was thought to be a short-scale phenomenon: locally parti- 
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cles get jammed and jailed in cages due to the interactions with the nearest 
neighboring particles in the liquid, resulting in a fast rattling motion but 
no structural relaxation. Stochastic nonlinear hydrodynamics is an effective 
theory for moderate and long length scales and therefore it was thought that 
many predictions were dependent on the short-scale cut-off or coincided with 
the previous one accidentally jS]. 

The second criticism is more recent and is related to the fluctuation- 
dissipation relation (FDR) between correlation and response function (see [T3| 
and next sections). The one-loop derivations of MCT presented in the litera- 
ture assume FDR in order to get the correct mode-coupling (MC) equations 
but actually they are incompatible with FDR. As a consequence they are not 
consistent and difficult to extend to more general cases. 

Despite these drawbacks the field theory approach is particularly appeal- 
ing and is gaining interest for many reasons. First, within this formalism it 
has been shown that MCT is a dynamic critical phenomenon and it leads 
to strong dynamical correlations in the four point density correlation func- 
tion ^H] (see for very important early insights based on the study of 
mean-field disordered systems). These results, which seem still out of reach 
of the projection operator formalism, have been indeed verified in simula- 
tions of model systems and quantitative MCT predictions fH] are in rather 
good agreement with numerical results [TZj(see also [IE])- Furthermore they 
provide an answer to the first criticism cited above: the MCT transition is 
related to growing correlations in the dynamics and is not a short scale phe- 
nomenon. The subtle point here is that the order parameter is the density 
correlator, a two point function, whereas in many other cases 1 the order pa- 
rameter is a one point function. Therefore the MCT equations have to be 
interpreted as mean-field equations on the order parameter. They play the 
same role as Weiss mean-field equations for ferromagnets. This makes clear 
that, first, they are very different from mode-coupling equations for critical 
phenomena. Indeed the latter are self-consistent equations of the correlations 
of the order parameter, and not of the order parameter itself. Second, to ob- 
tain the diverging correlation length one has to go beyond the mean-field 
MCT equations and compute the fluctuations of the order parameter, i.e. a 
four point density correlation function. From this perspective the fact that 
the order parameter, the density- density correlations, does not contain any 
diverging length is natural and does not imply that the MCT is a short-scale 

1 But not all! A famous example is the BCS superconducting transition. 
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phenomenon. 

The second reason which makes a field theory approach appealing is that 
it can be used naturally in off-equilibrium regimes. Indeed there is a strong 
interest in extending the MCT equations to regimes in which a liquid is aging 
or sheared. There have also been recently very interesting attempts to do 
that within the projection operator formalism [TTfl I2UI I2T] . 

Summary 

The aim of this paper is twofold: first we analyze field theories for the dy- 
namics of glass-forming liquids in terms of symmetries of the action. To 
our knowledge this has never been done and is very important in order to 
preserve physical symmetries in approximate treatments. In particular we 
will focus on time-reversal symmetry, which implies important relationships 
between observables. A well known example is the fluctuation dissipation re- 
lation (FDR) between correlation and response functions. Other important 
examples which have been discovered recently are the Jarzynski and Crooks 
equality for non-equilibrium processes [221 I2H] • We will focus on two types of 
field theories: the first is obtained from fluctuating nonlinear hydrodynam- 
ics (FNH) equations that describe in a coarse grained way the dynamics of 
liquids in which particles obey Newton equations j2U E] • The second one is 
obtained from the stochastic equations derived by Dean [22] • They are exact 
stochastic equations governing the evolution of the density field of interacting 
particles evolving with Langevin dynamics. The field transformations related 
to time-reversal symmetry that leave the action invariant are nonlinear in the 
fields. This is at the origin of the violation of FDR (and time-reversal sym- 
metry) by self-consistent perturbation theories that has been already noticed 
and discussed by Miyazaki and Reichman . We shall show that introduc- 
ing some auxiliary fields this problem can be solved and one can set up a 
self-consistent perturbation theory compatible with FDR. At this stage, we 
warn the reader that we will use the words "perturbation" and "perturba- 
tive" loosely throughout this article, as there will be no small parameter. In 
our context a "perturbation" is a formal series expansion, which radius of 
convergence may not be known. The second aim of this work is to analyze 
the field theory derivation of MCT using self-consistent one-loop approxi- 
mations which preserve the fluctuation dissipation relation between density 
correlations and response functions. This will clarify different issues related 
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to the field theory derivations of MCT and to extended MCT (and the re- 
lated cut-off of the transition). In particular we shall show that the two 
field theories assoctiated to FNH and Dean equations lead to essentially the 
same self-consistent one loop equations. These equations are very similar to 
the usual MCT ones but but with a slightly different wavevector-dependence 
of the vertex. This difference is clearly a drawback because the particular 
wavevector-dependence of the MCT vertex is very important for quantita- 
tive results. Furthermore, in the case of the Dean equation, our vertex leads 
to an apparently unphysical divergence. Our results show that some kind 
of resummation has to be performed in order to get MCT-like equations 
that are both quantitatively succesfull (as the one derived in 0) and respect 
explicitely time reversal symmetry. An interesting exact byproduct of our 
analysis concerns extended MCT. Indeed, we obtain the exact Schwinger- 
Dyson equation for the non-ergodic (Edwards-Anderson) parameter of the 
glass phase. We find that if the time-reversal symmetry is preserved by 
the perturbative self-consistent theory the MCT transition is not cut-off even 
when density and currents are coupled contrary to the conclusions of previous 
works (HI EE] • In those cases approximations violating time reversal symme- 
try produced a spurious cut-off of the dynamical transition. Our conclusion 
is that from a field theory perspective the mechanism which cuts the MCT 
transition off is a non-perturbative one unrelated to the presence of density- 
currents coupling and likely the same for Brownian or Newtonian dynamics 
as observed recently in numerical simulations |3D] and suggested also by very 
recent works jHSl HZj • 

This manuscript is organized as follows. In section I we introduce the two 
field theories for glass-forming liquids we focus on. From section II to section 
V, we focus on the field theory for the Dean equation, which provides a good 
illustration of the problems arriving with non-linear symmetries. In section 
II we show the fields transformations which leave the action invariant and are 
related to time-reversal symmetry. In section III we discuss how nonlinear 
field transformations affect perturbation theory and the origin of violation 
of FDR in the previous derivation of MCT. In section IV we show how one 
can restore FDR in self-consistent perturbation theories introducing more 
auxiliary fields. On the way an exact and compact structure of the Schwinger- 
Dyson equation is given, where the consequences of time-reversal symmetry 
have been exploited to the maximum. As an inllustration of the general 
strategy, we derive in section V mode-coupling equations which preserves 
FDR. In section VI, we present the result of the use of the strategy developped 
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in the previous sections to fluctuating nonlinear hydrodynamics. Section VII 
is devoted to a discussion of MCT related issues in particular the fact that 
the transition is not cut-off even in presence of currents. Finally, appendices 
deal with the most technical details. 

1 Dynamical field theory for glass- forming liq- 
uids 

The aim of this section is to introduce two field theories used to study glass- 
forming liquids such as supercooled liquids and dense colloidal systems. 

1.1 Brownian Dynamics and the Dean equation 

In the following we derive the field theory corresponding to a system of N 
interacting point particles evolving under Langevin dynamics which provides 
a good qualitative description of the dynamics of a dense colloidal suspension 
|26j . Actually one should also take into account hydrodynamic interactions 
due to the solvent but they will be neglected for simplicity. ^From a more 
theoretical point of view this is a limiting case in which, except for the 
density, there are no other conserved variables. We shall consider in the 
next section the field theory corresponding to Newtonian dynamics in which 
energy, momentum and density are conserved variables. 

The starting point is the Langevin equation which describes the dynamics 
of N interacting particles, evolving in Euclidean three dimensional space 
(coordinates are labelled by lower Roman letters, particles labelled by lower 
Greek letters): 

<9 t X a = - Vy(x a - X/3 ) + (a, (i) 

a</3 

where ( a is a Gaussian white noise, whose correlations are given by the 
Stokes-Einstein relation: 

(C a ,i(x, t)C/3j(x', 0) = ZTSapSij 8(x - x') 6(t - f). (2) 

Here, V is the pair potential between particles and the time is expressed in 
units of the microscopic diffusion constant. 
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The potential V is denned up to an additive constant. We shall tune this 
constant so that f c/ 3 x V(x) = ~~ which will be convenient later 2 . Observ- 
ables measured in experiments are functions of the local density p(x, t) = 
^ a 5(x — x a (t)). Using Ito calculus, Dean has shown that the density of a 
system of particles obeying (JTJ) obeys the following Langevin equation 



<V(x, t) = V • Lo(x, t) v ^7y) + >7( x > *)> ( 3 ) 

where 77 is a Gaussian random noise, whose correlators are: 

Mx, tMx', 0) = 2Tp(x, t)V • V'5(x - x')5(t - 0- (4) 

The prime on the gradient means that it acts on functions of x'. The oc- 
currence of multiplicative noise here is not surprising, as the local density 
does not fluctuate in empty regions. The density functional T for a fluid of 
average density p can be written as: 

T\p\ = T j rf 3 x p(x) (]n^--l\+^J rf 3 x J d 3 x' p(x)p(x')V r (x - x') 
= -TS[p\+T int [p\- 

(5) 

The Langevin equation (jHJ) is equivalent to the equation derived heuristically 
by Kawasaki j27j, when — /3V(x) is replaced by the direct pair correlation 
function c(x) , and the functional defined in (jSJ) becomes the Ramakrishnan- 
Youssouff (RY) density functional 

The dynamic average of an observable *4 over thermal histories can be 
expressed as a functional integral over the density field: 

{A) = J VpA[p](5[d t p^,t) - 77] (x, *)]>„, (6) 

where R[p, 77] (x, t) is the RHS of ©, (•}?, means the average over the Gaussian 
noise 77 and 6[-] is a functional delta function. This is the standard procedure, 
see jSHj, to derive field theories from stochastic equations. The only subtlety 
is the absence of a Jacobian. This is due to the fact that stochastic equations 

2 More precisely, we put the system in a box of volume f2, replace the potential V(x) 
by V(x) - ^ (j^ - / d 3 xK(x)j and take the limit fi -> 00. 
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with multiplicative noise (J3J) are defined following the Ito prescription and 
therefore the Jacobian is a constant that can be absorbed in the definition 
of the functional integral. This is also related to the Markov property of the 
stochastic differential equation © (in the Ito discretization). 

By using an integral representation of the functional Dirac distribution 
through a conjugated field p and averaging over the noise 77, the Martin- 
Siggia-Rose-de Dominicis-Janssen (MSRJD) action is obtained p?T] : 

(A) = JvpJ VpA[p] e?M, (7) 

with 



S[p,p] = Jd 3 *J dt |p(x,t) 

+ Tp(x,t)(V/3(x,t)) 2 }, 



p(M)V-* 



flp(x,t) 



(8) 



or explicitly: 

S[p,p] = J d3 *J dt <p(x,t) 



4p(x,t) + TVV(x,t) (9) 



+ V ■ (p(x, t) j d 3 x' VV(x - x')p(x', t^j 



Tp(x,t)(Vp(x,t)) 5 



For clarity, and as the dynamical action involves one-time quantities, we 
shall not write the explicit time dependence of the fields in the rest of the 
paper, except for correlation functions involving fields at different times. For 
instance p(x, t) will be written p x , and J <i 3 x J dt will be replaced by J . 

1.2 Fluctuating Nonlinear Hydrodynamics 

We shall now recall the field theory used to investigate the dynamics of 
compressible liquids in particular close to the glass transition. In this case 
particles evolve under Newtonian dynamics and the derivation of the corre- 
sponding stochastic equations is not from first principles contrary to Dean's 
equation. The resulting equations are meant to be a generalization of hy- 
drodynamic equations to intermediate time and length scales, and are not 
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expected to lead to an accurate description of the physics occuring on short 
time and length scales. Hence, they describe the evolution of slow variables 
(associated to conserved quantities) subject to a thermal noise (correspond- 
ing to the fast degrees of freedom that have been integrated out) and they 
have been called Fluctuating Nonlinear Hydrodynamics equations (FNH). A 
phenomenological derivation and a discussion of the equations can be found 
in 0123]. In the following we will focus on the FNH equations used by Das 
and Mazenko jH] to investigate the problem of the glass transition. They fo- 
cused only on density and momentum as conserved variables but in principle 
the energy can be introduced as well. The equations of Das and Mazenko 
read: 



where p and ^ are respectively the density field and the ith component of 
the momentum density field and rji is a white Gaussian noise with variance 
(r/j(x, t)r/j(x', t')) = 2TY i j(x.— x!)5(t— t'). The effective free-energy functional 
is T = Tkin + Fu, with 



where po is the density of the system, m is the particle mass and c(x) is the 
direct correlation function. Note that the potential term, Tu, coincides with 
the Ramakrishnan-Youssouff functional. The Poisson brackets {■, ■} and the 
Yij are chosen so that the continuity equation for the density is verified and 
in such a way that the linearized equations coincide with the usual linear 




(10) 



(11) 




(12) 
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hydrodynamic equations: 
r fi (x-y) 



-Vo ( 3 ViVj + (JyV 



CoViV,- 



<J(x 



(13) 



= L#£(x - y) 
{p(x),^(x')} = -V^(x-x')p(x) 
{^(x),p(x')} = -p(x)V^(x-x') 

{^(x),^(x')} = -Vi<5(x - x'Mx') -^(x)V^(x-x'), 

and ?7o and Co are respectively the bare shear and bulk viscosity Using the 
previous definition one can rewrite the equations (|1(JI12|) in a more explicit 
way: 



dtPx 
dt9i, x 



-ArV 



(14) 



5p2 



u 



9i,x9j, 
Px 



: )-X>fe)+** < 15 > 



Following the strategy of the previous sections we find that the corresponding 
MSRJD action is the integral of: 



~Px 

9i,x 



d t px + Vi p x 



<5pz 



+ Tgi :X Lijgj 
' $9j,x 



(16) 



V7 ^ 

9j,xViT 1" 

ogj, x og jtX 



where p and g are respectively the conjugated field used to express the delta 
functions corresponding to Eqs. ()14I15|) . and summation over indices is im- 
plicit. In the next sections, we will focus on BDD only. We shall show later 
in section El how to generalize the results to the case of FNH. 



2 Time-reversal symmetry and fluctuation dis- 
sipation relations 

Time-reversal symmetry relates the probabilities of a path and of its time- 
reversal counterpart in configuration space. This is a very important symme- 
try obeyed by systems in thermodynamic equilibrium and it has far reaching 
consequences. In particular all physical correlation functions are invariant 
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under time-reversal, correlation and response function are related by the 
fluctuation-dissipation relation. In the context of dynamic field theory time- 
reversal symmetry is related to a transformation of the fields leaving the 
action invariant. In the following, we first recall the standard field trans- 
formation for one particle evolving under Langevin dynamics. We will later 
show how this can be generalized to the more complex field theories intro- 
duced in the previous section. 



2.1 Langevin dynamics 

Let us now focus on one particle evolving under Langevin dynamics with 
additive noise. Denoting the particle position X(t) and the external potential 
V the Langevin equation reads: 

d t X(t) = -VV(X(t)) + V (t), (17) 

where rj is a Gaussian white noise with zero mean and variance 2T. After 
the introduction of a conjugated field X, the dynamical 
action reads 



SM}= / dt 



oo 



CO 



-X(t)-(d t X(t) + VV(X(t))) + TX(t) 2 , (18) 



where the dynamical field is <p = (X,X). It is easy to check that So is 
invariant under time-reversal: 

° : { x(t) ->x(t)-±d t x(t). (19) 

As a consequence correlation functions of X(t) are invariant under this fields 
transformation and therefore they are invariant under time-reversal. The 
application of the fields transformation to the response function, which can 
be shown to be equal to (X(t)X(t')) , leads to 

(X(t')X(t)} = (X(t)X(t>)) - ^d t (X(t)X(t')}. (20) 

For t > t', the response vanishes, and we get 

(X(t')X(t)} = ~d t {X(t)X(t')) } (21) 

which is the the FDR between correlation and response functions. 
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2.2 Brownian dynamics for the density field 

We will now focus on the stochastic evolution of the local density field for 
interacting Langevin particles (in what follows we use the acronym BDD). 
Note that the extension of the ideas in the next paragraphs to generic mul- 
tiplicative noise is straightforward. First we derive the expression of the 
response function in terms of fields average and after we show two fields 
transformations which lead to a physical time-reversal symmetry and, hence, 
to FDR. 

2.2.1 Response function 

In order to establish FDR, one has to write the expression for the response of 
the system to an external potential. Here, as we probe density fluctuations, 
this force is an external potential, which is taken into account by adding 

an extra term J- ex t[p] = — PxPx to the free energy (jHJ. The response 

J X 

Rxx' (t, t') at time t and position x to an infinitesimal external force switched 
on at time t' and position x' is defined by 

(p(x,% = (p(x,0) M=0 + J d 3 x" J dt" R xx ,(t,t")fi(x",t") + o(/i), (22) 

where is the average taken with the free energy functional T + T ex t- 
Thus, expanding the path integral to first order in p, one gets: 

(p(x,t)) At = (p(x,t)) At=0 

+ J d 3 x" I™ dt" (p(x, t)p(x", t") V • (p(x", t") Vp(x", t")) . (23) 

However, due to causality the time t" in the integral runs until t only, and 
integrating twice by parts, one obtains the expression of the response: 

Rxx>(t,t') = -<p(x,t)V- (p(x',t')Vp(x', *'))>• (24) 
This response function is the same as the one studied in 



2.2.2 First expression of time-reversal symmetry 

The first field transformation related to time-reversal symmetry is 

T : ( J ~* 7* - (25) 

I Px ->• Px + fx, 
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where / verifies 

V • (p x Vf x ) = ~d tPx . (26) 

The field / plays a role similar to the longitudinal part of a current for the 
density field. By integrating twice by parts, the variation of © under T is 
found to be 

~J d t F\p x \. (27) 

Hence, at equilibrium the action © is invariant under T. The density field 
is affected in a simple way by T: p(x, t) — > p(x, —t). This implies that any 
average of the type p( x «) ti)) is invariant under time-reversal. By making 
the change T in the expression (|24|) of the response function and using the 
condition (f2Hj). we get: 

7jPt'C xx i{t — t') = R xx '(t — t') — R xx '(t' — t), (28) 

where C xx >(t — t') = (<5p(x, t)5p(x', t')). This is the FDR, which can be easily 
generalized to correlators of more density fields: 

^dfC Xl ... XnX i (t\ , t n: t ) Rx\---x n x' (^1 ' i t-nt t ) (29) 

Rx\---x n x' ( ^1 ) tny t ), 

where 

Rx 1 -x n x'(h--- ,t n ,t') = -(V- (p(x',t')Vp(x / ,t'))n'°( X "^)) ( 3 °) 

1=1 

and 

n 

C xv .. XnX '{t 1 • ■ ■ ,t n ,t') = (p(x',t') J^p(xi,ti)). (31) 

i=l 

2.2.3 Second expression of time-reversal symmetry 

As said above, the time-reversal symmetry may be expressed through another 
transformation of the fields. Indeed, consider the following change: 

op x 
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We remark that this transformation involves only the non-physical field p, as 
does T. Thus, both transforms clearly give rise to the same relations between 
correlation functions. 

As above the dynamical action can be easily proved to be invariant under 
U. The time-reversal symmetry of any density correlation functions follows 
directly. FDR can also be derived using the identity: 

Wx - i) ^ > = °- < 33) 

which implies 

(p(x, t)<9,p(x', t')) =(p(x, t) V • (p(x', t')Vj^-^j ) 

-2T(p(x,t)V.(p(x',t')Vp(x',t')))- 

Splitting the last term in two identical parts and applying the transformation 
hi to one of them leads us to 

<p(x,t)«9,p(x',t')> =<p(x,t)V- (p( x '^')V^^y)) 

-T(p(x,t)V-(p(x',OVp(x',t'))) 

+ T(p(x, -t)V ■ (p(x', -t) Vp(x\ -f ))) 1 J 

Mp(x,-t)V.(p(^-OV^^)). 

The two terms containing no p in the RHS cancel and one gets (J28|) . 

^From U one can derive another useful identity as follows. Applying the 
transformation to G xx >(t — t') = (p(x, t)p(x', t')) one gets: 

(36) 

As we shall show in paragraph 14. 61 this identity allows substantial simplifica- 
tions when T is quadratic. Indeed in that case U is a linear transformation 
and (}36j) then becomes a linear relation between the density-density corre- 
lator {6p(x,t)5p(x',t')) and the naive response Q xx >it — t'), which are both 
two-points functions. 
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3 Nonlinear symmetries and perturbation the- 
ories 



As will be emphasized all along this article, a physical symmetry is related 
to fields transformations that leave the action invariant. In the following we 
shall show from a general point of view why the linearity of these transforma- 
tions is essential to construct symmetry preserving perturbative expansions. 
None of BDD and FNH field theories has this property, namely the fields 
transformations related to time-reversal symmetry are nonlinear. In the fol- 
lowing we shall highlight all the complications which arise in perturbation 
theories because of the nonlinearity, focusing on BDD. In particular this will 
make clear why mode-coupling approximations (MCA) generically violate 
time-reversal symmetry and hence FDR. 



3.1 General discussion 

For the sake of generality, let us consider a field <j) a (where a = (x, t) in our 
case) and a generic field theory: 



n 



cxp 



(37) 



where g is the coupling constant controlling the perturbative expansion, V 
is the interaction and the field <fi is assumed to have a zero mean. The bare 
propagator is G and the bare action Sq[4>] = — |0 q (Gq ) a ,/3<f>p + gV(4>). 

Now consider a linear transformation O of the field (fi a — > <f)' a = O a ^^ 
which leaves So invariant. We require that the integration measure is left 
invariant (the contrary giving rise to so-called anomalies in quantum field 
theory), which means that O must have determinant of modulus equal to 
one (or at least to a constant). Then we have the following Ward-Takahashi 
(WT) identities: 

(0ai =O ai ,/ 3l ---O an7 p n {(f)p 1 ---(f) l 3 n ). (38) 

One immediately sees the problem when the symmetry is not linear. Consider 
for instance a transformation of the form <fi a — > cf)' a = O a ^<^. The WT 
identities are 

(0ai-"-0a„) =O ai j 3l -'-O anj i3 n (<l?p 1 '--<f?p n ). (39) 
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Clearly, a nonlinear symmetry induces relations between correlation functions 
of different orders in the fields. A concrete example is the proof of FDR 
obtained in the previous sections for BDD field theory. The general idea 
behind mode-coupling approximation and other self-consistent approaches 
is to provide closed equations for two point functions. Hence, because the 
transformation mixes correlation functions of different order in the fields it 
seems difficult, if not impossible, to construct MC-like approximations which 
preserve the symmetry. In order to understand where perturbation theories 
fail in the case of non linear symmetries it is useful to recall what are the key 
ingredients that make them work when the symmetry is linear. 

Let us first focus on standard (non self-consistent) perturbation theory. 
A linear symmetry implies a relation between two point functions, see Eq. 
()38[) . Why is this relation preserved order by order in perturbation theory in 
gl The reason is that the potential part V(<p) is itself invariant under O and 
thus flHHj) is also true when both LHS and RHS are computed at any order in 
g. For concreteness, let us focus on the MSRJD field theory for one particle 
evolving under Langevin dynamics. The dynamical action can be split into 
two parts: Squad + Sint, with 

/oo 
dt -X(t)-d t X(t)+TX(t) 2 (40) 
-oo 

and 

/oo 
dtX{t)-VV{X{t)). (41) 
-oo 

Both parts of the action are invariant under time-reversal and one can expand 
in powers of g: 

oo / x n 

(X(t')X(t)) = J2^-(X(t')X(t)S INT [<!>} n )GAus S , (42) 

where (■) gauss stands for the Gaussian average with action Squad- However 
Sint is itself invariant under the action of O, thus the identity (|21jl . FDR, 
is true at any order in perturbation theory. 

Let us now focus on self-consistent perturbation theory. This approxima- 
tion scheme consists in cutting the perturbation series at a given order and 
replacing the bare propagator by the dressed one in the formal expression 
of the self-energy (some diagrams present in bare perturbation theory have 
to be neglected in order to avoid double counting). Doing so one obtains 
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self-consistent equations for the evolution of the dressed propagator which in 
principle can be improved considering more diagrams. Note also that these 
self-consistent equations can be often shown to be exact in some large N 
limit, where the field has become a N component field 

The formal way to obtain these self-consistent perturbation theory is 
through the Legendre transform T(G) = log F(E(G),G) of two points func- 
tions which generates all the two particles irreducible diagrams |32| : 



exp T(G) = f Y[d(paexp - - ^ <l>a(G l ) a ft<l>p + gV(4>) 

a I a,/3 



a — G/3 } a) 



a;3 



(43) 



where £ is determined by the condition ||j = 0, i.e. £ such that the propa- 
gator of the theory equals G. De Dominicis and Martin have shown that V 
is equal to 

- l -Tr (Gq 1 G) + ^Tr InG + <f> 2PI (G), (44) 

where 4>2Pi{G) is the sum of all 2PI diagrams with the full propagator G 
used as internal lines [221 • The physical propagator of the theory is obtained 
by finding the G which makes T stationary, i.e. solving f§ = 0. Hence, it 
is immediate to check that the self-energy as a function of the propagator 
is S(G) = 2 ^qq 1 ■ Practically the usual self-consistent or mode-coupling 
approximations keep the lowest non-trivial diagrams, the bubbles, in <p2Pi ■ 

Let us now prove that the symmetry is preserved if the field transforma- 
tion c() a — > (p' a = O ai p(j)p that leaves the action invariant is linear. For the same 
reason as before we assume | det 0\ = 1. Under this transformation the prop- 
agator transforms into G' a p = M a> p ;a i t p>G ' a t t pi where M a ^ a i t pi = O a ^ a iOp y p< 
(|detM| = 1). We shall prove that this relation is preserved in the self- 
consistent perturbative expansion. 

The first thing we want to prove is that T(MG) = T(G) which in a certain 
sense is expected since a physical symmetry should not change the value of 
the functional. This is easy to prove starting from the expression for T(MG) 
and rewriting (p a = O a> p(j)'p and changing variable in the functional integral. 
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The measure as well as the action is invariant so that one finds 



-iZapM&^h + gVte) (45) 

Calling £' = SM and noticing that dr (^ IG "> = is equivalent to l9r ^,' G ^ = 
we obtain that indeed T(MG) = T(G). Furthermore since this is true for any 
value of g, the coupling constant, then this identity is true for T expanded 
to any given order in g. The consequences of this identity are particularly 
useful. First if G is a solution of = so is MG. Thus, if there is a 
unique solution - the symmetry is unbroken - then G = MG. In the case of 
time-reversal symmetry this identity leads to FDR. The fact that the solution 
is unique is expected in our case since time-reversal symmetry is certainly 
not broken at equilibrium. Since this is true for any value of g then it is 
true when T is expanded to any finite order in g, i.e. for the self-consistent 
equations written to any finite order in g. Finally the last remark following 
from T(MG) = T(G) is the symmetry relation verified by the self energy: 
S(MG) = S(G)M. 

To recap, we have shown that a symmetry related to a linear transfor- 
mation in the field is a symmetry of the self-consistent equations written to 
any finite order in g. That is the reason why FDR is preserved for standard 
MSRJD Langevin field theory. Note that a very different proof can be found 
in . 



F(MG) =log J JJ# a exp 



3.2 Subtleties in perturbation theory 

In this section we focus on BDD field theory and we highlight the difficulties 
and the failures of perturbation theories vis a vis FDR. The case of FNH 
is conceptually identical but practically more clumsy because of the larger 
number of fields. Let us start with bare perturbation theory. In order to do 
a perturbative analysis it is convenient to separate the Gaussian zero mean 
part of the local density field from the interacting one by introducing density 
fluctuations in Q Sp x = p x — po- This gives: 

S = / (s 0lX + s INTiX ) , (46) 



21 



with: 



s , x = p x (-d t 6p x + TV 2 5 Px + p jf V{x - *')V 2 5 Px )j + Tp (Vp x ) 2 

stnt,, = T8p x (Vp x ) 2 + p x V ■ Up x £vy(x-x')^ . (47) 

This may be written in a more compact form through the bidimensional 
vector field (Sp,^. In Fourier space, the inverse of the propagator of this 
field is: 



r-i = ( iuj + Tk 2 (i+p Po v(k)) \ 

~ V -iu} + Tk 2 {l + {3 Po V{k)) -2T Po k 2 J 



(48) 



However, including the potential in the form (J48|) makes it practically impos- 
sible to write weak coupling expansions that preserve time-reversal symmetry. 
Hence we shall drop it out of the inverse of the propagator: 



! / C (k,u) G (k,u) \ 1 
\G* (k,u) ) 



(49) 



iu + Tk 2 

-iu + Tk 2 -2T Po k 2 , 



and we shall treat the quadratic term as an insertion in a line. This gives 
the following Feynman rules: 

• bare density correlator: 
Co(k,uO 

• bare naive response: 
« Qo(k,u) 

• line insertion: 
-*o*- - Po k 2 V(k) 

• potential vertex: 

X k" I (k • k'V{k>) + k • k'V(k")) = T(k, k', k") 
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• noise vertex: 




-Tk' ■ k' 



The bare density correlator is 



C (k,cu) 



2T Po k 2 



(50) 



u 2 + (Tk 2 ) 



and the bare naive response 



£o(k,cj) 



1 



(51) 



Tk 2 — iu 



Due to the form of the vertices, diagrams with tadpoles and hence corrections 
to the average density vanish to all orders (the momentum at the entrance 
into the tadpole is zero). 

In terms of density fluctuations, the response is 



R xx >(t,t') = -(5p(x,t)V- (5p(x',t')Vp(x',t'))> -p V 2 (<5p(x,t)p(x',O) 



where Q is the naive response and x is an "anomalous" response. Having a 
look at (|52|). one sees that part of the non-triviality of the FDR arises from 
the anomalous response x> which itself comes from the multiplicative aspect 
of the noise, or equivalently from the nonlinearity of the transformation of 
the fields associated with time-reversal. 

A perturbative expansion in powers of the potential preserves the sym- 
metry. Indeed the dynamical average of any function A[p, p] may be written 
as 



where Spree contains the terms of the action which do not contain the 
potential. The key point is that both Spree and Sy are invariant under T, 



Xxx'{t,t') - po V*G xx >{t,t'), 



(52) 




(53) 
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and thus FDR is verified to any finite order in the expansion in power series 
of Sy. However Spree is not quadratic in the fields. The nonlinearity of the 
fields transformation relate the noise cubic term to quadratic terms. Thus 
this vertex must be taken into account non-perturbatively. In other words, 
to compute the correlators at a given order in the power series expansion 
in Sy, one has to include contributions at all orders in the noise vertex. 
This is a difficult but not impossible task because, at a given order p in Sy, 
the diagrammatic expansion contains a finite number of diagrams, due to 
the absence of the propagator connecting two p's. Indeed, if one considers 
a diagram with p potential vertices and q noise vertices contributing to a 
correlation function of r p's and s p's, one must have q + s < p + r. Hence 
such a diagram must have less than p + r — s noise vertices 

Thus bare perturbation theory which preserves the nonlinear symmetry 
can be set up but is considerably more complicated than the usual one. 
Finally, this discussion makes it clear that any approximation which drops 
a part of the diagrams of the full expansion in terms of the noise vertex is 
expected to be in contradiction with the FDR. This is indeed what happens 
in self-consistent approximations as we shall show below. 

3.3 Violation of fluctuation dissipation relations in self- 
consistent perturbation theory 

Let us now focus on self-consistent perturbation theory in particular on the 
mode-coupling approximation introduced by Kawasaki |33] that consists in 
neglecting vertex renormalization. First, we write the Schwinger- Dyson (SD) 
equations 

G Q 1 ■ G = 1 + E ■ G, (54) 

where £ is the self energy: 

^ <»> 

and the associative product ■ is defined as follows: 

POO 

(A-B)(k,t) = dt' A(k,t-t')B(k,t'). (56) 
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Causality and reality of the density auto-correlator imply that the self- 
energies verify: 



£ P p(k,t) = £^(k,-t) (57) 
= 0. (58) 

The first diagrams contributing to the self-energies are 



,(2) 



^(2) 






+ 





Diagrams of higher orders all contain vertex renormalization. Hence, if one 
neglects renormalization of both vertices, the SD equations (J54)) become MC 
equations for Q (for t > 0): 



d t g(k,t) = 
d t C(k,t) -- 

with: 

£pp(k, t) 



- Po Tk 2 (1 + (3 V(k)) G(k, t) + dt> £ pp (k, t - t') C(k, t') 
- Po Tk 2 (1 + /3 V(k)) C(k, t) + J Q °° dt> E p/5 (k, t - f) Q(k, t') 

+ J^ oo dt'^ p (k,t-t')C(k,t'), (59) 



d 3 q 
(2^)3 
d 3 q 



^>pp(k, t) 



(2tt) 3 
d 3 q 



(2tt) 3 
d 3 q 



(2tt) 



£(k, t) C(k - q, t)r(q, k, k - q)r(k, q, k - q) 

G(k, t) g{k - q, t)q • (k - q)r(q, k, k - q) 
C{k,t) C{k-q,t) 2 (60) 
Ug(k, t) C(k - q, t)q ■ kr(q, k, k - q). 



These equations are not compatible with FDR, as can be seen from the solu- 
tion of SD equations at low orders in the potential. FDR is trivially verified at 
order zero. At order one the MC equations are exact and hence they are auto- 
matically compatible with FDR. Incompatibilities appear at order two, where 
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diagrams such as those shown in Fig.Q the first diagrams contributing to ver- 
tex renormalization, have to be taken into account in the non self-consistent 
perturbation theory and therefore also in the self-consistent one (as discussed 
previously, in order to preserve FDR one has always to take into account the 
contributions at all orders in the noise vertex contribute to FDR). This sug- 
gests that if one wants to improve the approximation by keeping for instance 
the first vertex corrections, one has to include at least all the diagrams of 
order two (in powers of the potential) in the self-energy. In that case, the in- 
compatibility with FDR would be an effect of order three. However, nothing 
guaranties that the violation of FDR by the self-consistent approximation 
is attenuated when the order of the approximation is increased. Another 
consequence of practical importance is violation of causality in (|59|) . where 
times in the integrals are not restricted in [0, t}. On the contrary, when time- 
reversal - and thus FDR - is preserved, times in the integrals run from to 
t. To conclude this paragraph, we remark that incompatibilities with FDR 



Figure 1: Example of diagram which contribute at order V 2 to FDR but is 
absent of the MC equations. 

have arisen from an explicit breaking of the time-reversal symmetry which 
is due to the nonlinearity of the field transformation related to time-reversal 
symmetry. 

4 Restoration of time-reversal symmetry in 
perturbative expansions 

In this section, by introducing some extra fields, we construct a generalization 
of the field theory described previously in which time-reversal symmetry 
corresponds to a field transformation which is linear. As a consequence all 
the problems described in the previous section are eliminated. In particular 
this allows to set up self-consistent perturbative equations which preserve 
FDR. We explain in detail our procedure for the BDD field theory ; the 
result of this strategy for the FNH field theory will be given in section El 
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4.1 Introducing extra fields 

As already discussed in order to overcome the violation of time-reversal sym- 
metry in perturbative expansions we need to make the related fields trans- 
formation linear. This can be achieved introducing some auxiliary fields. 
Furthermore we also introduce a response field such that the response func- 
tion becomes explicitly a two-fields correlation function. We have found two 
different fields transformation related to the same physical symmetry (time- 
reversal). As a consequence we can make each one or both of the symmetries 
linear which leads to three different field theories. In this section we shall 
treat the case in which only either the U or T transformation is linearized. 
The case of the completely linearized theory (both symmetries are made 
linear) is considered in Appendix E. 

Let's consider the time-reversal transformation T first: we start from the 
identity: 

(A) = J VpAlp](5(d t p{x,t)-R[p,r)](x,t))), (61) 

where R[p, rf\ (x, t) is the RHS of 0. We now plug into the functional integral 
the representation of the identity 

/ v f n 6 ( v ■ G°* v -k) + Tfbp*) ( det [- y • ^* v )]) = l ( 62 ) 

We put a minus sign in such a way that the operator inside the determinant 
is positive definite. Thus we don't need to take the absolute value of the 
determinant. We exponentiate the delta function using an auxiliary field / 
and we introduce fermionic fields (f> and to exponentiate the determinant. 
As a consequence there is a new term to add to the previous action that 
reads: 

J fx fv^V/A + ^d tPx \ - J p x V(f> x ■ V4> x . (63) 

Furthermore we introduce also the field ip = V ■ (pVp) which allows for the 
usual two-fields correlator representation of the response function. This leads 
to introduce a conjugated field ip for the Fourier representation of the delta 
function related to i/j. The final action is the integral of 

S x = - PxdtPx + IpxT Tf) x 1p x + lp x (lp x - V, • (px^xPx)) 

-/I ^ N _ («) 

+ fx ( j^dtPx + V x ■ (p x V x f x ) j - p x V x ■ ((j) x V x (f) x ) . 
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This action now remains invariant up to the boundary terms under the fol- 
lowing linear transformation 7^: first invert the time t — > —t then change the 
fields in sequential order as follows: 



Px * Px fx 

Ipx -> Ipx + \d t px 

ipx -> 4>x + Tf x (65) 

fx -> -fx + Tf x + ^j x + Tp x 

<. fx * /a; 



Formally we write this as <ft — O ■ <f), where <fi = (p, p, ip, ip, f, f, <ft, <j)) T and 
<j>{x.,t) = 0(— x, t). This implies an identity for correlators: 

G = O ■ G ■ T . (66) 

The transformation has a determinant of modulus one, as a product of simple 
transformations with this property. 

Let us now show that this transformation implies FDR. Consider R xx >(t' ' — 
t) = — (p(x, — t)ip(x.' , —£')). Under the transformation 7{ this transforms into 
Rxx'it — f) — \dt'C xx i{t — t'). Thus the equality (jSSj) implies in particular 
Rxx'(t — t') = R xx /(t' — t) + j;dt>C xx i(t—t') which is the fluctuation-dissipation 
relation. 

We now show how to linearize the second transformation U. We intro- 
duce the field 9 = 5J r /5p and the conjugate one 9 to exponentiate the delta 
function 3 . 

The action S is then transformed into an integral of 
s x = -Pxdtp x + Ti>AY,P,) 2 + X (Ox - ^p) - Px(V x p x )(V x 9 x ). (67) 
The associated linear transformation U\ is 

U J Px > ~Px + ^p9 x (68) 
\ 9 X 9 X — 7fd t p x , 

As before we write it as <p = O ■ 0, where <fi = (p, p, 9, 9) T and 0(x, t) = 
</>(— x, t). The same identity ()66|) holds for correlators. Again the transfor- 
mation has a determinant of modulus one. The response function can be 

3 The usefulness of introducing these two fields when dealing with the BDD field theory 
was noticed by C. Chamon and L.F. Cugliandolo from a sligthly different perspective |28j . 
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written as R xx >(t — t') = (p(x, t)#(x', t')). Using the transformation U\ we 
find again that (p(x, — t)#(x', — £')) equals R xx >(t — t') — ^d t 'C xx '(t — t'), hence 
FDR. Note that this second field theory is considerably simpler than the 
previous one because it has less fields. 

One might wonder how the final result depends on the choice of the 
linearized field theory. As far as self-consistent perturbation theory and 
MCT is concerned we have written the dynamical equations obtained by 

• a) using the completely linearized theory where both symmetries are 
rendered linear. 

• b) only the fields involved in the transformation T\ in addition to p. 

• c) only the fields involved in the transformation IA\ in addition to p. 

• d) writing the terms of the action such that the potential V is in one 
of the vertices (in the goal of making the link with standard MCT). 

At the order of one loop, we have found the same sets of equations for corre- 
lation and response function at long times in all cases. This is not surprising, 
since the different transformations do not affect the physical fields and change 
the response fields in the same way. We thus expect this to be valid at all 
orders. Indeed, as we will show below, one gets closed equations for the 
dynamical evolutions of correlators involving only the fields p and 9. In ad- 
dition, this tends to confirm that FDR makes the results robust with respect 
to the choice of extra dynamical variables. Thus in the following, we will 
focus on the simplest theory written above in terms of p, p, 9 and 9 only, the 
choice of the fields to work with being merely a matter of taste. We refer the 
reader interested in the theory with all fields introduced above to Appendix 
E. 

4.2 Minimal theory preserving fluctuation dissipation 
relations and its basic properties 

In order to avoid cumbersome calculations, we shall from now on describe a 
minimal (in the sense of the number of fields) theory for which the symmetry 
associated to FDR is linear. It is the theory produced by U\. In the following 
we shall describe the WT relations for correlation functions and self-energies 
due to time-reversal symmetry. These are particularly useful because they 
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make clear that there are only three independent correlation functions or self- 
energies out of sixteen. In the next section we shall write down the general 
form of the Schwinger-Dyson equations. Using the WT relations we find that 
these reduces only to three independent integro-differential equations. 

Henceforth we use the density of action (|57jl and the symmetry U\. We 
split the density of action into a Gaussian part s an d the interaction part 

SlNT'- 

so, x = -pxd t 5px - Tp p x V 2 x px + X 9 X - 9 X (W * 5p) x + p a p x V 2 9 x {m) 



Si NT,-. 



T5 Px {Vpxf - S Px (V x p x )(V x 9 x ) +T9 X J2- 



n>l 



Spx 

Po J 



where 



and 



(W*p) x = [ W(x-y)p y 

J y 



Whc) = V(x) + —She). 

Po 



(70) 



(71) 



(72) 



It is useful to remark here that each term of the expansion in n is indepen- 
dently invariant under U\. 

One gets the following relations for correlators from identity (|66|) applied 
to the case of tl\ (the limit t, t' — > oo, r = t — t' is taken): 



Cpp^xx' 

C p e,xx' ( r ) 



T 



AT 



9(r) 



Cee,xx'{T) 



T 

ll-{Q{r)C^ xxl {r)) 
9(r) d 



rp Cpp,xx' \T) Rxx' 1 



T 



(73) 
(74) 
(75) 
(76) 



with the obvious notation C a b jX x'(t — f) = (a x (t)b x /(t')). These identities 
imply only three independent two-fields correlators, namely C pp>xx r,C p Q !XX i 
and Cee,xx'- Moreover, the causality of the theory, which is insured by the 
Ito discretization, is explicit. A perturbative proof of causality can be found 



30 



in Appendix B and comes from the causality of the bare propagator, which 
in Fourier space is: 



/ 



2T P0 k 2 



1 



(w 2 +pgfc 4 vy(k) 2 ) 
1 

-iu+p k 2 W(k) 

2T P0 k 2 W(k.) 
(w 2 +pgfc 4 W(k) 2 ) 
Pok 2 

\ -iu>+p k 2 WQf.) 



iu+p k 2 W(\i) 



W(k) 
' iuj+p k 2 W(k) 





2Tp k 2 W(k) 
(u; 2 +p 2 k 4 W(k) 2 ) 

W(k) 

-iw+p k 2 W(k) 
2Tpgk 2 W(k) 2 
{u) 2 +p 2 k±W{\t) 2 ) 



Pok 2 
iuj+pok 2 W(k) 



-kj 

iu>+p k 2 W{k.) 





\ 



/ 



(77) 



iu)+pok 2 W(k.) 

The bare propagator helps also to understand the anomaly in (|73j) where the 
time derivative of the Heaviside function is present. One can easily see that 
the bare propagator (|77jl has the form: 

C 0$(k, uj) = — 1 + function of(k, uj). (78) 

It is easy to prove using (J69J7()|) that there are no diagrammatic corrections 
to the constant part of the above correlator, that is in the above equation 
the constant part persists perturbatively and we can write: 

9(r) d 



C' 



(79) 



One can find (see appendix C) that C p g tXX >(0) = T. The latter together with 
(J73J) allows to write (f7H|) . 

Recall that the Schwinger- Dyson equations have the form Gq -G = 1+D, 
with D — E • G. Thus, S transforms under IA\ in the following way: 

£ = 0~ T ■ £ • C^ 1 , (80) 

enforcing the following constraints on self-energies: 



%(k,r) 

Spe(k, r) 
S^(k,r) 
The other elements 



^£fc(k,r) 



6(r) d 

T 

9(r 



E ee(k, r) 



(81) 
(82) 



£ PP (k,r) (83) 

k,-r)-re(r)S^(k,r). (84) 

and See) vanish. One can write the equation 



T 

TQ(-r)E 



for Ep P (k, r) in the following form: 



^pp(k, 



l__d_ 

fdr 
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e(r)E^(k,r) 



(85) 



As before we start with the diagrammatic analysis recovering a set of dia- 
grams which contribute to the 5-function term. Their direct resummation 
is cumbersome. We used the Schwinger-Dyson equations to identify the 8- 
function term. Note that as for correlation functions there are only three in- 
dependent self-energy terms. All the others are either zero or related through 
WT identities. 

4.3 Dynamical equations 

In the following section we transcribe the full dynamical equations and leave 
their derivation for Appendix C. There are only three independent equations. 
Using the transformation laws under IA\ it is easy to see that there are not 
more than four independent equations. The proof that one of this equations 
is trivially verified once the other three are verified is more tricky and done 
in Appendix D. 

Let us first choose r > and write the equations which will give the time 
evolution of the correlators at strictly positive time- difference. The values at 
t = of the correlators will be obtained from the study of the singularities 
(remember there is a 5(t) in the RHS of (|54|)). 

We first consider (G^-G — E-G) ( g p (k, r) = 0. The corresponding equation 

is 



d T C pp (k,r) + p k 2 C pe (k,r) = / dtZ P 8(k,T-t)C pe (k,t) (86) 



Now consider (G 1 • G — E ■ G) p g(k, r) = 0. The corresponding equation is 



Finally consider (G • G — E • G)§ (k, r) = 0. The corresponding equation is 




d T C pe (k, r) + p k 2 C e e(k, r) = -E^(k, r)C>(k, 0) 
+ / dtZpe(k,T-t)C e g(k,t)+ I dtZ §e (k,T-t)d t C pe (k,t). 



(87) 




+ 
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As discussed at the beginning of this section, there is an extra equation that 
at first sight might seem independent from the first three (but in fact is not, 
see Appendix D). It comes from (G^ 1 •(? — £• G)§ e (k, r) = and reads: 

W(k)C pd (k, t) - C w (k, r) = Is w (k, 0)C>(k, r) - i^(k, r)C pd (k, 0) 

+ / dtX §e (k,T-t)C B e(k,t)-^ f dtZ §§ (k,T-t)d t C p e(k,t). (89) 
Jo 1 Jo 

We stress here that obtaining such consistent structure - which is the 

exact one - for the dynamical equations is far from trivial. Approximations 

which violate FDR generally lead to a different structure. In such cases, 

the resulting conclusions about the existence of the glass transition are very 

suspicious because they might be due to the violation of this fundamental 

structure. 



4.4 Static limit 

The static equations are obtained by taking r = in the SD equations. This 
is an advantage of the present field theory approach, in which the statics is 
included in the dynamics, compared to derivations of dynamical equations 
based on Mori-Zwanzig formalism. We remark that the closed set of static 
equations cannot be obtained from the evolution equations as written above, 
as the singularities at r = have been excluded from the latter. The correct 
derivation is given in Appendix C. The correlators and their derivatives at 
initial time are: 

C>(k,0) = T (90) 
C pp (k,0 + ) = - Po k 2 C>(k,0) = -Tpok 2 (91) 

C ee (k,0) = W(k)C pd (k,0) =TW(k) (92) 
C>(k,0 + ) = -^(k)Tp k 2 -T£^(k,0) (93) 

CUM) = ^ + _^E w (k J 0)C7 /v (k ) 0) > (94) 
where C stands for d T C. The equation ()94)1 is obtained from ()88)) and ()107)) . 



4.5 Equation for the non-ergodicity parameter 

In MCT, the glass phase is characterized by a non-zero value of the so-called 
non-ergodicity parameter |3T| I38j . which signals the existence of infinite-time 
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correlations. The non-ergodicity parameter /(k) is defined as follows: 



In MCT, /(k) is zero in liquid phase and jumps to a non-zero value at the 
glass transition. Here the equation for the non-ergodicity parameter can be 
obtained by taking the limit r — > oo in (j86H89|) . 

We first give an argument in support of the vanishing of C p g(k, oo) and 
Cgg(\t, oo). The physical interpretation is the following. We label by a each 
ergodic component into which the system may break (if no ergodic to non- 
ergodic transition occurs then there is only one such component). Then in 
each component the system decorrelates completely at long times, and we 
can write: 



where each component is characterized by its weight W a and the static av- 
erage (-)q, inside it. However, the average force (^) Q vanishes, and hence 
C p g(k, oo) too. This also applies to Cgg(k,oo). We stress here that the 
possibility to express infinite time averages in term of static quantities is 
intimately related to the fact that the dynamical equations are consistent 
with the dynamical symetries at equilibrium. Indeed this guarranties that 
the asymptotic measure is Gibbsian. Let us add that C p g(k, oo) = may also 
be seen as a direct consequence of C pp (k, oo) = 0. Furthermore, Eq. 
gives Sg e (k, oo) = 0, and putting (|9*3Jl into (J9*3j) leads to the exact equation: 



Of course this equation is too difficult to be analysed directly because E^ 
contains an infinite number of terms that cannot be resummed. The same 
occurs in the Mori-Zwanzig formalism, where E^ is related to the long time 
limit of the memory kernel. One has to resort to approximation schemes, 
the simplest one being replace E^ by its one loop expression. That will be 
discussed in the following. 

4.6 Quadratic density functionals 

In this section we discuss the particular case in which the density functional 
is taken to be quadratic in the fields. This corresponds to make a simple 




(95) 




(96) 




(97) 
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quadratic expansion of the entropic part S and to replace in the FNH case 
the kinetic term by g 2 /2p where po is the average density. The former 
case has been already studied in ^3] and the latter in [SH]. In these works, 
especially in [T2|, the issue of preserving FDR has been already investigated 
and it has been shown that one-loop mode-coupling equations preserve FDR 
in this case (contrary to the general case). 

^From our perspective the reason of this preservation is simple: in these 
cases the fields transformations become linear even without introducing extra 
fields. Let us consider in detail the BDD case studied by Miyazaki and 
Reichman [TB*] . 

When expanded at order two, the entropic part is written as 

S[p] » S[p ] -T f M. (98) 

Jx Z P0 

Thus the free energy functional J-\p\ becomes quadratic: 

F\P\ = \I PxPx'W (99) 

* Jx,x' 

The equilibrium measure from (|9*9*j) is Gaussian: 

If d 3 k p(k)p(-k) 

v[p]oce -2j Wf^W (ioQ) 

S(k) = (5p(k)5p(-k)} 



WW 

where we have used J c? 3 x W(x) = 0. This Gaussian form considerably sim- 
plifies the structure of the resulting field theory. Indeed, the transformation 
14 defined earlier becomes linear: 

U ' { I - -Px + (3f x , W(* - x') Px , (101) 

As stated previously, the use of U makes it possible to derive an identity 
between the density- density correlators and the naive response. One can 
write 

(p(x, t)/3(x', t'))=- (p(x, -t)p(x', -t')) 

(102) 



+ (3 J d 3 x' W(x-x') (p(x,-t)p(x',-t / )), 
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or equivalently in Fourier space: 



g(k,Lu) + g(k,-u) 



C(k,u) 



(103) 



The symmetry IA is linear and thus does not mix correlators of n fields with 
correlators of n + 1 fields, unlike T. This makes the identity (jl()3|) valid at 
any order in perturbation - self-consistent or not - as explained in section I. 
In Appendix A, in order to give a concrete view of how diagrams have to be 
put together to lead to ()103J) . a diagrammatic perturbative proof if given, as 
well as the derivation of a similar relation between self energies: 



This identity makes ()103|) compatible with the SD equations Using 
these identities, it is possible to show that FDR is preserved by self-consistent 
one-loop theory and also to all order in self-consistent perturbation theory 
as explained in Appendix A. 

A final remark about the quadratic case is that from (J97|) . it follows 
f(k) = 0, as there is no vertex involving 9 and hence must vanish on 
general grounds. 

5 Mode- coupling approximation preserving fluc- 
tuation dissipation relations 

Now, we will carry out the MCA for the BDD field theory that we studied 
in the previous sections. The corresponding analysis of fluctuating nonlinear 
hydrodynamics follows the same guidelines (see section^. It is important 
to stress the difference between our and previous approaches. We set up 
a self-consistent diagrammatic expansion that preserves automatically time 
reversal symmetry and FDR. As a consequence the corresponding MCA equa- 
tions, or any other approximated expression of the self-energies, will do that 
too. 

5.1 Mode-coupling approximation 

We first focus on the term J x x ^j^-. The entropic part of T gives a contri- 
bution J x 6 X log ( 1 + — J to this term. In order to compute the dynamical 



Epp(k, t) 



1 



E pp (k,t), t>0. 



(104) 



S(k) 
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partition function one may expand the logarithm in powers of 2fe This gives 

an infinite number of vertices of the type T<=p. (£)'. A crucial point 

is that, as stressed above, all the terms in the action coming from different 
powers in the series expansion are independently invariant under the trans- 
formation U\. We also can put some couplings in front of the other vertices 
and carry out truncated expansions at different orders for different vertices. 
However as we focus on the one-loop theory, and for simplicity, we shall treat 
these vertices as if they were all of order T. Then there are two ways of deal- 
ing with the vertices arising from the expansion of the logarithm. On one 
hand one can take into account all these vertices. However this leads to the 
sum of an infinite number of terms, whose meaning is not clear, due to the 
presence of an infinite number of tadpoles which contribute to static vertex 
renormalization. On the other hand one can truncate the series expansion of 
the logarithm. However one needs to go beyond first order in order to take 
into account nonlinearities. For simplicity we cut the series at order two, but 
the calculation can be in principle extended to any order. 



5.2 Expression of the self-energies 

Within the approximation described in the previous paragraph, the self- 
energies read: 

/3, 



S <^ (M = J ^)^pl ) Cpe ^ t)Cpp ^~^ t] (105) 
S^(k,t) = ±y , ^|(k.q) 2 C' w (q,t)C' /¥ ,(k-q I t) (106) 



-(k-q)[k-(k-q)]C p0 (q,t)C pe (k-q,t) 
<i 3 q 



S wOM) = j^J j^C pp {q,t)C pp {k-q,t). 



(107) 



It is instructive at this stage to see how the MC equations derived by Miyazaki 
and Reichman ^2] in the case of the quadratic expansion can be obtained 
from the above equations. The simplifications arising from the linearity of U\ 
will then become clear. Using the method used to derive (J102)) . one obtains 

W(k)C pe (k,T)=TC ee (k,r), (108) 
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and (|1U2|) can be simply written as 



W(k) C pp (k,r)=TQ(r) C>(k, r+) + T6(-r) C>(k, (-r)"). (109) 

Using these two identities one can eliminate C P and Cog in (|SH|). In addition, 
when the entropy is expanded at quadratic order, all self-energies except Epg 
vanish, and the equation which remains is the one derived and discussed 
in [T3|. 

5.3 Static limit 

The static equations are obtained by taking r = in Eqs. ()105H107ft . Re- 
markably, it is identical with the one which would be obtained by doing the 
same MCA on the static theory involving the density functional T[p\. Indeed, 
after truncation of the entropy at cubic order, the density auto-correlator is 



JV 5pp(y)p(y')e-^ 
JVSpe-PKlSp] 



Cpp,yy'{t - Uj r^rmz^r^ , l nu J 



with 



f 3 [5p} = iy rf 3 x^^ - pj d?*j d 3 x'lU(x - x')5p(x)5p(x'), (111) 

where we have used the fact that J <i 3 x W(x) =0. It can be easily verified 
that using the MCA as described above for the correlator defined by (jllOj) 
gives again (|9*^jl. 

5.4 Equation for the non-ergodicity parameter 

Using the one-loop expression ()107|) in (}9~Tj) . one gets the equation for the 
non-ergodicity parameter: 

/(k) 1 f d 3 



1 " /(k) 
Writing this in the form 



^ / 7^V(q)/( k " q)5(q)5(k - q)S(k). (112) 
2 PoJ ( 27r ) 



/(k) = l- — -, (113) 
m(k) 
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one sees that /(k) = and /(k) = 1 are two solutions of (|112|) . However 
the second solution, which in addition makes the integral in m(k) diverge, 
implies C pp (k, t) = 1 for all t and thus has to be rejected. We stress that in 
contrast with standard MCT the g-dependence of the vertex (in particular, 
the absence of the term (k • qc(q) + k ■ (k — q)c(k — q)) 2 ) does not inforce the 
convergence of the integral in (|112|) . However, the equation (I94j) which gives 
the structure factor is also ill-defined as the integral over q diverges too. We 
have solved numerically these equations putting a cut-off for large values of 
|k|. We have found a |k| dependence of /(k) very similar to the usual one of 
MCT. However, as the cut-off goes to infinity our numerical solution seems 
to converge toward the solution /(k) = 1 for all k. However, although there 
is a clear physical cut-off to the description in terms of Langevin equations, 
the cut-off dependence we found is clearly unphysical. It might be that 
convergence is only obtained via a further resummation of diagrams that 
renormalize the vertex. We leave this problem for a future work. Ignoring 
this cut-off problem (which is in a sense solved in the FNH case) for the time 
being, our MC equations have the following properties: 

• as in the standard MCT, one can derive a schematic theory [37| by 
assuming that the structure factor is dominated by a single mode 
S(k) = 1 + AJ(|k| — fco). Then the schematic equation for the non- 
ergodicity parameter becomes identical to that of the schematic MCT. 

• We have also checked that the dynamical critical properties are the 
same as for standard MCT. 

• The static equations can be reduced to an equation involving the den- 
sity field only, which is identical to that we can get from the equilibrium 
density functional by doing equivalent approximations. The way the 
theory has been written thus makes it possible to treat the dynamics 
in a similar way to the statics. 

6 Fluctuating Nonlinear Hydrodynamics 

In this section, we describe the case of FNH. As the derivations follow closely 
those of BDD, we just give the results. 
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6.1 Time-reversal symmetry and fluctuation-dissipation 
relations 

This time there are four response functions produced by the extra term T ext = 

~ J x (PxPx + gx ' Px) 



(p(x,i)) M =(p(x,t)) AJ=0 + J d 3 y J dt"R pp>xx ,{t-t")n{y,t") 
+ J Sy J dfR k pgtXX ,(t-t") Pk (y,t") + o(v,p) 



(114) 



<<fc(x, *)>„=o + / d 3 y I dt"R% )Xa ,{t-t") V i{y,u 



+ d 3 y dt"Rl PjXxf (t-t")fi(y,t") + o(ii,p) 



which gives: 

Rpp,xx' (t ) t ) 



Rgp,xx' {t 

7>Sr 



(115) 



Rpq,xx' if t ) 



Rqq,xx'(t t ) 



(p(x,t)V-(pg)(x',i')> (H6) 
<<fc(x,t)V-(pg)(x',t')> (H7) 
(p(x,t)p(x',t , )V fc p(x',t')) + (p(x,t)^(x',t')V fc &(x',t')> 
+ (p(x,t)V i (^)(x',t')) " (p(x,t)L fc ^(x',t')> (118) 
{0 fc (x,i)p(x,t)V,p(x',t')) + (y fc (x,^ i (x' J OV, &(x , ,0> 
+<< /fc (x,t)V i (^)(x , ,t , )> " (^fc(x,t)L K ^(x , ,f)> (119) 



The transformation of the fields associated to time-reversal is now: 



V : ^ 



Px 



-t 



1 5F 



(120) 



Spa 
1 bJF_ 
T8 Sx 



This transformations leaves the action invariant up to boundary terms: 



dig;, + —a t p x + Vj — 7 — 

og x dp x \ dg itX dg jtX 



(121) 
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Following the procedure used for the BDD field theory one can see that 
density correlation functions are invariant under time-reversal and also derive 
FDR. We stress again that the naive self-consistent perturbation theory for 
the FNH violates time-reversal symmetry in the same way it does in case 
of the BDD. An extension of the original model is needed providing linear 
time-reversal symmetry in order to satisfy the symmetry in perturbation 
expansion. 



6.2 Restoration of time-reversal symmetry in pertur- 
bative expansions 

In order to make the transformation V linear we introduce the two additional 
fields 9 — t£ and v = We are lead to add 

dp <5g 



^ 

OPx 



- v a 

.1- 



8T 



- / Ox 

J X 

to the action, which becomes the integral of 

-px [dtPx + Vj (p x Vi)] + Tg itX L i:j g jtX 

!!,.,- [d t 9i,x + Px^iO + Vj (gi, x Vj) + //,., V;r, + LijVj] 

8F~ 



Vx — T 

dpx. 



The corresponding linear time-reversal transformation is: 



t - 


-> -t 


§x 


~* Sx 


Vx - 


-+ "Vx 


Px - 


-> -Px 


§x 


§x 


@x - 


-+e x + 


Vx " 


-> "Vx 



\dtPx 
- \d&x 



(122) 



(123) 



(124) 



As before, the action is conveniently split into Gaussian and interaction 
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parts: we expand all the terms in powers of Sp/po. 



S 0,: 



-p x (d t Sp x + p V ■ Vj,.) - g i)X [d t g ijX + p Vi8 x + LyVj^] (125) 



-l'!li.,-I-ij>lj..r - 9 X 9 X + 6 X (W FNH -k 5p) x - w x ■ v x H v x • g x 

Po 



SlNT, X 



-p x V ■ (5p x v x 

9i,x9j,x^i^j,x 



n>0 



- Sp x 



- g itX Vj(g itX v 



'ST Tin M 

"x p 

p m (f Q 



*)5p x + n6 x -^- 



Sp 



n-l 



Po +1 ' 



with 



W 



FNH, 



X = — 



T 

m 



-*(x) 
Po 



c X 



?72 



(126) 



(127) 



This expansion in powers produces two series of vertices, each being inde- 
pendently invariant (also order by order) under the field transformation V±. 

Different equalities between correlators (and self-energies) result as con- 
sequences of the use of Vi as for the BDD field theory. In particular using 
that the response function is the correlation between p and 9 one gets FDR 
between correlation and response. The set of all relations and the simplified 
dynamical equations are presented in Appendix F. 

6.3 Static limit 

The analysis of singularities of the SD equations at short time difference 
gives: 



C>(k,0) 
C gv (k,0) 

W(k)C pp (k,0) 
C ee (k,0) 



Po 



C gg (k,0) 



C vv {k, 0) 



T 
T 

T + is^(k,0)C pp (k,0) 
TW(k) 

T+^E^(k,0)C gg (k,0) 

T 
Po' 



(128) 
(129) 

(130) 

(131) 

(132) 

(133) 
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6.4 Equation for the non-ergodicity parameter 



We now focus on the equation for the non-ergodicity parameter. As for BDD, 
due to long time decorrelation inside ergodic components, only C pp , C pg ± and 
Cg± g ± (where g -1 is the transverse current) can have a non-zero limit when 
t — > oo. This is confirmed by the analysis of the dynamical equations given 
in Appendix F. We however do not expect frozen currents, and thus it is 
quite reasonable to asume that C P9 (k, oo) and C gg (k, oo) do vanish. Since at 
least one of those appear in any diagram of S^(k, oo), E^k, oo), S-^(k, oo), 
Sg{)(k, oo) and S^~(k, oo) those self-energies also vanish. As a consequence 
one obtains the non-perturbative equation for the non-ergodicity parameter 
using (|234|) and its static limit (jl3U|) : 

--U(k)X^(k,oo). (134) 



1 - /(k) T2 

This structure is identical to the one found for BDD and is an exact result. 
Any general approximation (one loop, two loops, etc) for the self-energy 
on the right hand side, lead to a non linear equation on /(k). As we will 
discuss in detail later previous works have obtained very different structures 
because in those cases time-reversal symmetry was violated. This may be 
very dangerous because it can generate spurious results as it is indeed the 
case for the cut-off of the transition. For example in the analysis of [Hj this 
modifies strongly the general structure of the Schwinger- Dyson equations 
and implies that the non-ergodic parameter has to vanish. However, this has 
nothing to do with the physical mechanism that cut-off the MCT transition 
and is just an artefact of having violated the time reversal symmetry. 



6.5 Mode-coupling approximation 

We restrict ourselves in this section to an approximation similar to the one 
used for BDD, truncating the vertex series in 5p/po at the lowest order: 

sint,x = -pVV • (5p x v x ) - 5p x g x ■ V9 X - gi, x Vj(gi, x Vj, x ) 

-gi,xgj,xViVj, x - v x • gx^r - x7 -^ - t^—^x^t- ( 135 ) 

Po 2 Po 2m Po 

We do not write the set of all equations at the one-loop level. We just remark 
that as in the BDD case the static correlation functions can be obtained from 
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the dynamical equations. They coincide with the ones obtained doing the 
analoguous approximation on the static theory. At the level of one-loop, the 
expression of H§§ reduces to the one previously obtained for BDD (up to a 
multiplication by the mass). We want to stress that it is a coincidence which 
is absent when higher orders schemes are considered. We believe that this is 
due to the fact that the static effective free energies are different in the two 
field theories. 

As a consequence the MCA equation for the non-ergodic parameter is the 
same as the one obtained for BDD and the previous remarks applied also to 
this case. We have also checked that the critical long-time behavior of the 
correlation function close to the MCT transition is the standard MCT one 
and coincides with the one obtained in the BDD case. Note that contrary to 
BDD there is now a cutoff to regularize the integral over q because FNH is 
not valid on short-length scale. 

7 Relation with previous works and issues re- 
lated to Mode- coupling Theory 

In this section we would like to put our work in the context of field the- 
ory derivations of MCT and discuss what can be learned from the non- 
perturbative structure of the equations which we derived, and from the re- 
sulting mode-coupling equations. 

The first issue that we want to discuss is the field theory derivation of 
mode-coupling equations. One can find in the literature [SHI E] different field 
theories derivations of the original full fc-dependent mode-coupling equations 
introduced and studied by Gotze et al. |H E] • All these derivations are not 
consistent because they assume some identities (related to time-reversal sym- 
metry) that are incompatible with the same self-consistent equations used to 
derive MCT. Indeed the Kawasaki derivation of MCT using field theory [SH] 
starts from a BDD field theory (in which the potential is replaced by a term 
proportional to the static direct correlation function). Kawasaki computed 
for the original (without extra fields) BDD theory and assumed that E^ p 
was related to it by a symmetry transformation similar to FDR for correlation 
functions: 

S^(r) = -ia^r). (136) 
Looking at the whole set of one-loop equations it is easy to see that this rela- 
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tion is not verified. Actually the structure of the Schwinger-Dyson equations 
is such that if this relation was verified there would have been a fluctuation 
dissipation relation between the correlation function and C pp (t) which is not 
the response function as discussed previously. 

We now focus on the Das and Mazenko jH] field theory derivation of MCT. 
In the original papers and all subsequent ones there is a linear relationship 
that has been assumed to hold between C pp and C pp (see eq. 6.62 of 0). This 
relation, that has been used to close the equations on C pp , is true only in the 
hydrodynamic limit (long time and length scales) or for a purely quadratic 
free energy density functional. Instead, time-reversal symmetry generically 
implies a more complicated identity that is the generalization of (|36|) to the 
FNH case: 

*)|tt^> = ©(* - (p( x - *)p(x', 0> + - *) <p« Op(x, t)>, 

dp(x', V) 

(137) 

where T is the free energy functional introduced for the FNH field theory. 
Forcing the relation between C pp and C pp is useful to close the equations but 
(1) it is very dangerous as already discussed in 16.41 and as it will be explained 
later and (2) it is inconsistent with time reversal symmetry that imposes a 
different relation, eq. (|1H7|) . 

Contrary to these two approaches our derivation of MC equations pre- 
serves time-reversal symmetry because it is constructed upon a self-consistent 
expansion that preserves automatically this symmetry. This is important for 
the self-consistency of the theory and it is crucial in order to study off- 
equilibrium dynamics as discussed in the introduction. At one loop, it leads 
to a vertex with a different fc-dependence than usual MCT. This is problem- 
atic for two reasons: first the good quantitative results of MCT depend cru- 
cially on the vertex fc-dependence and, second, our vertex leads to a strange 
and cut-off dependent behavior or the mode-coupling equations. Our conclu- 
sion is that more refined (e.g. higher loops) self-consistent approximations 
have to explored is such a way to obtain MCA that preserves time reversal 
symmetry and is at the same time quantitatice succesfull as the one devel- 
oped in jH]. The investigation of this issue is left for future studies. 

An interesting consequence of our derivation and our equations concerns 
extended MCT and the cut-off of the transition. Using field theory jS] and 
also projection operator formalism jU] equations beyond standard MCT have 
been obtained. These are called extended MCT equations and have been con- 
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jectured to contain the "hopping or activated" effects that destroy the MCT 
transition. A striking result is that these "hopping or activated" effects are 
present only if there is a density- current coupling. As a direct consequence, 
they would be expected for particles evolving under Newtonian dynamics 
(and thus described by FNH) and not for particles evolving under Brownian 
dynamics (and thus described by BDD). Recent numerical simulations have 
shown [ID] that quite the contrary happens: the same effects cutting off the 
transition are present for Brownian or Newtonian particles. It has been also 
believed for some time that ideal MCT transition might occur (i.e. is not 
avoided) for Brownian dynamics but it is now clear that it is impossible since 
there are mathematical theorems proving jJTj that the self-diffusion coeffi- 
cient of hard spheres evolving under Brownian dynamics never vanishes at 
any finite temperature or chemical potential. 

As a conclusion the Mode Coupling Transition is expected to be always 
replaced by a cross-over. The crucial question is what is the physical mech- 
anism that cut-off the transition, how capture it in an analytical theory and 
whether it is the same or not for Newtonian and Brownian dynamics. In 
our derivation we find at long-times the same transition, and at one-loop the 
same equations, in the Newtonian and Brownian cases, hence there is no sign 
that the cut-off of the transition is due to the presence of coupling to currents. 
This conclusion remains valid at any finite order in the self consistent per- 
turbation theory (one loop, two loops, etc.). In practice any approximation 
will provide an expression for S^(k, oo) as a function of /(k) that has to be 
plugged into eq. ()134j) or its BDD counterpart. There is no general argument 
implying that the resulting non-linear equation on /(k) cannot have a non 
vanishing solution and indeed one expects that the transition present at one 
loop carries on at higher loops. That has been indeed verified in toy models 
|4*2] as Langevin particle inside a double well potential. Our conclusion is 
therefore that the cut-off mechanism is due to non-perturbative effects that 
cannot be captured at any finite order of the self-consistent perturbation 
theory. 

This is in clear contradiction with previous field theoretical works that 
found a cut-off mechanism for one loop self-consistent expansion. In the 
following we unveil, using exact results, that that mechanism was only due 
to the approximations that were used and that violated explicitely violation 
of time-reversal symmetry. 

The evidence for an avoided transition in the field theory derivation comes 
mainly from two works. Das and Mazenko |U] found that when all one- 
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loop diagrams are considered the transition is cut-off. However, this was 
due to the relation they assumed between the two correlation functions C pp 
and Cpp. In fact, as we have explained earlier, C pp cannot have a plateau. 
Therefore by forcing this relation one kills artificially any possibility of having 
a glass transition. Another way to put it, as discussed in 16.41 is that this 
relationship alters completely the non-perturbative structure of Schwinger- 
Dyson equations. Whereas it is not possible to conclude anything about the 
cut-off of the transition just looking at the general equation for the non- 
ergodic parameter, see section 16.41 the relationship assumed in jU] between 
Cpp and Cpp plus the general form of the Schwinger- Dyson equations imply 
that the non-ergodic parameter has to be zero. Indeed, consider one of the 
exact Schwinger- Dyson equations derived in the Appendix F: 

d T C pp (k, t) - ip kC vp (k, r) = — [ dtS^(k, r - t)d t C pp (k, t) 

1 Jo 

+ f J Q dtE p z(k,T-t)d t C gp (k,t) - ^ dtZ pP (k,T-t)C 9p {k,t) 

- 7fj dtZ pS (k,T-t)C vp (k,t) (138) 

and the exact fluctuation dissipation relation derived in Appendix F: 

C^(k,r) = ^C>(k,T). (139) 

The linear relationship assumed in j2] between C pp and C pp implies that 
if Cpp has an infinite plateau at the transition so does C pp and, using the 
FDR relation ()139j) . C p g. However, a plateau of C p g is incompatible wit the 
exact equation (|138|) . Indeed, integrating eq. (|138|) over r between zero 
and infinity, we would get an infinite right hand side and a finite right hand 
side 4 . Therefore one would conclude that the first hypothesis, i.e. an infinite 
plateau for C pp , has to be wrong. However, this is only due to having forced a 
relation between C pp and C P p\ Thus, the cut-off of the transition found in [§] 
is spurious. The relation between C pp and C pp is valid only on hydrodynamic 
length and timescales. 



4 Only the third term on the right hand side would give an infinite contribution since 
C vp cannot have a plateau. This can be found by inspection studying the other equations 
or just on physical grounds because at the glass transition amorphous circulating currents 
are not expected. 
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We finally remark that other arguments complementary to ours have been 
given recently [22] > in favour of rejection of the cut-off mechanism derived 
within extended of MCT provided in 

The other work that is cited as a supporting evidence for the cut-off of 
the transition is the one of Schmitz et al. jSHl- However, they considered a 
purely quadratic free energy functional. We have seen in our derivation that 
the transition comes from the density nonlinearity term in the free-energy 
functional and therefore it is absent in the quadratic case. In [36J FDR is 
respected and the corresponding equations are a particular case of the ones 
we derived. The fact that there is no transition is the natural consequence 
that theories with quadratic free energy functional are too simple to lead 
to any MCT transition in the FNH as well as in the BDD case. In fact, 
contrary to the case of non quadratic functionals, there is no transition at 
any finite order in the self-consistent perturbation theory simply because the 
self-energy £^ in eq. ()134|) is identically zero! Note that in the Brownian 
dynamic case this has been already found, although in a different way, by 
Miyazaki and Reichman [T3] . 

Conclusion 

In this paper we analyzed field theories for the dynamics of glass forming 
liquids focusing on their symmetry properties. In particular we have shown 
that straightforward perturbation theories generically do not preserve time- 
reversal symmetry. We have found that time-reversal symmetry is related 
to a nonlinear field transformation that leaves the action invariant. The 
nonlinearity is the source of the problems. Introducing some auxiliary fields 
we have shown how to to set up perturbation theories that preserve time- 
reversal symmetry and hence FDR. 

Furthermore we have critically revisited the field theory derivations of 
MCT showing that they always assumed some relations which are actually 
incompatible with the self-consistent equations. Our derivation is completely 
consistent and preserves FDR but leads to a different vertex than in the 
usual MCT. This leads to similar qualitative results but clearly different 
quantitative results. Furthermore strange and spurious divergences appear 
if we do not put an ultraviolet cutoff on momenta integration. We leave for 
future work an accurate investigation of this puzzle that will certainly need 
the introduction of some kind of diagram resummation. 
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We have also reconsidered the evidence for the cut-off of the MCT transi- 
tion when a coupling between density and current is present. We have shown 
using exact results that there is no obvious cut-off mechanism which acts or- 
der by order in self-consistent perturbation theory. ^From this perspective 
there is no difference between Brownian and Newtonian dynamics. Whether 
or not density-current couplings are present, we have found the same formal 
structure for the equations on the non-ergodic parameter in the glass phase. 
Actually at one loop the corresponding equations for the non ergodic param- 
eter are identical for BDD and FNH. This structure is fundamentaly different 
from the one previously obtained [5] that suggested a cut-off mechanism only 
for Newtonian dynamics; the latter was due to the assumption of a linear 
relationship between correlation functions, that although valid in the hydro- 
dynamic limit, is not verified in general and that, via the general structure 
of the Schwinger-Dyson equations, force the non-ergodic parameter to be 
zero. The correct relation, eq. (j!37|) . has no influence on the existence of the 
MCT transition. We conclude that the MCT transition has to be cut-off by 
non perturbative mechanisms. A very recent work |37j shows by a schematic 
approximation that once any factorization approximation is avoided, i.e. at 
a non perturbative level, the MCT transition indeed disappears. 

Despite the problems - due to excessive simplicity of approximations such 
as the MCA considered here - related to our vertex we think that our ap- 
proach is promising since it automatically connects statics to dynamics in a 
precise way. It remains the issue of finding appropriate approximations for 
the self energies. However, one has not to worry without about compati- 
bility with time-reversal symmetry, which is guarranted by the form (j86H89J) 
of the Schwinger-Dyson equations. This step was missing in the previous 
attemps to derive MCT equations from field theory, where approximations 
were done too early. Thus, it opens the way to dynamical equivalents of the 
self-consistent schemes developed in liquid theory |J3J|2]. This was already 
initiated in a pioneering work $Jd\ where the relationship between dynamic 
and static theory for the glass transition was investigated. 

Finally we remark that the identification of field transformations associ- 
ated to time-reversal symmetry may help developing other non-perturbative 
approaches. For instance in the non-perturbative renormalization group lUIj 
approach, an Ansatz for the effective dynamical action has to be made. The 
choice among different possible Ansatze is very wide but can be drastically 
reduced by the use of symmetries. We hope that the analysis of the symme- 
tries performed in this work may be helpful in applying NPRG techniques to 
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the the glass transition problem. 

Acknowledgements We are deeply indebted to J. -P. Bouchaud for very 
useful discussions, a careful reading of our manuscript and his participation 
to the early stage of this work. We also thank L. Berthier, C. Chamon, P. Cal- 
abrese, L.F. Cugliandolo, K. Miyazaki and D. R. Reichman for discussions, 
in particular K. Miyazaki and D. R. Reichman for a careful reading of this 
manuscript. A. L. acknowledges support from European training program 
STIPCO under grant HPRN-CT-2002-00319. AA and GB are partially sup- 
ported by the European Community's Human Potential Program contracts 
HPRN-CT-2002-00307 (DYGLAGEMEM). 

Appendix A. Perturbative proofs for the quadratic 
theory 

A.l. Proof of (TTUSft 

In this appendix a perturbative proof of (|104J) will be given (in the quadratic 
expansion of the free-energy). The identity (|1(J3|) will follow from the relation 
between self-energy and connected two-points functions: 

G~ l - G 1 = E. (140) 

Here the order n in perturbation theory consists of diagrams with n vertices 
B[p,p], the vertex including two parts B = B\ + B2, B\ and B2 being the 
noise and W vertices respectively. 

In this appendix, we will always have t > 0. As already said above, the 
splitting of p into po an d 5p leads to the splitting of both noise and interaction 
vertices in two parts in the dynamical action. The first one, linear in p is 
quadratic and can be seen as a "mass" insertion, whereas the second one is 
cubic and gives the vertices used in the perturbative expansions. However it 
turns out that both parts (the one with po and the other) are independently 
invariant under the symmetry U. Thus there are two strategies do carry 
out perturbative expansions. One consists in keeping the quadratic term 
proportional to po in the bare propagator, while the other consists of treating 
it as a mass insertion. For the purpose of this appendix, both strategies are 
almost equivalent, and thus we will follow the first one which generates less 
diagrams. In this case, the bare propagator is given by PSJ) . At the bare 
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level (|1U3|) is verified, that is 

£o(M) + £ (k, -t)=pW(k) <7 (M). 



(141) 



Causality implies that for t > t' the terms of order n + 2 (n > 0) in the self 
energies may be written formally as S = Q±pi, with 

eg^k, * - f ) = J 0/ *(p + p' - k) *(q + q' - k) 
xr(p, k, p') <( 4r(q, k, q') (Sp(p, t)6p(p', t)B[5p, p} n p(q, t')Sp(q', t>)) 



-2Tq ■ q'(5p(p, t)Sp(p', t)B[Sp, p] n p(q, t')p(q>, *')> 



(142) 



9g +2) (k, t-t') = J^J 0/ 0/ *(p + p' - k) S(q + q' - k) 
xr(p, k, p') J 2r(k, q, q') (8p(p, t)8p(p>, t)B[8p, P ] n 8p(q, t')5p(q>, f)) 



-4Tk • q(Sp(p, t)Sp(p', t)B[Sp, p} n p(q, t')5p(q>, t')> } . 



(143) 



The subscript 1PI indicates that only one-particle irreducible diagrams are 
considered. In term of diagrams 







pp 




e 



pp 




where the ZVs are obtained by doing all possible contractions in the averages 
of ()142j) and ()143|) . The set of diagrams with m vertices will be noted QS m \ 
The exponent ^ will in general design classes of diagram of order m exactly, 
while the subscript ( m ) will include all orders until m. It is convenient here to 
consider the subset of 1PI diagrams of order m as the subset of all diagrams 
of order m, minus the classes of connected but not 1PI diagrams (subset 
^4( m )) and of disconnected diagrams (subset B^). We will then proceed as 
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follow. First we will prove that the components of verify an identity similar 
to (|1(J4|) . and that the same holds when restricted to class B^ n \ Then we 
will prove simultaneously that (jlU4|) is true at order m and also holds when 
restricted to class A^ m \ The sum of all diagrams belonging to class is 
noted A (m) . 

Applying U in all terms with p explicitly written in (|142j) and (|143|) . one 

gets 

e£ +2) (k, t - t) = J faf ^ r S(p + p' - k) % + q' - k) 
xr(p, k, p') (5p(p, -t)5p(p', -t)B[5p, p] n 5p(q, -t')5p(q', -t')) 

^ q .w(k)+ q . q ^(qO w ^ _ 2/?q . q/ w ^ W(q')) (144) 

e£ +2) (k, t - t') = / faf *frf 8(p + p' - k) % + q' - k) 

xr(p,k,p') (2 k - qI ^ (q)+ 2 kqW(q,) (Sp(p,t)5p(p',t)B[5p,prSp(q,t')5p(q',t')) 
-4k • q W(q)(5p(p, -t)Sp(p', -t)B[5p,p} n 5p{q, -t')5p(q', -t'))) . (145) 

All correlation functions involving p(q, —t') or p(q', —t') have vanished due 
to causality Indeed they contribute to P p(k, i! — t) — 0p P (k, t — t') which 
vanishes for t > t'. Finally, using the fact that correlators involving only 
B[5p,p] n and p are invariant under time-reversal, the identity ()104|) is found 
at order m = n + 2. 

The same steps can be follows straightforwardly for disconnected dia- 
grams, as it amounts to splitting B[5p,p] n into several powers of B[Sp,p] 
which contribute to different connected components. The crucial point here 
is the invariance of each term B[Sp,p] under U. 

Now we prove that diagrams of order n verify the identity ()1U4|) holds 
at order n and also within class A^ n \ To do so, we proceed by induction. 
Assuming that and verify (|1U4|) for all p < n — 2, we will show that 
£^ and A^ verify IflMjl . 

The starting point is however n = 2. At this order: 

9 (2) (k,t-t') = S (2) (k,t-t') (146) 
A {2 \k,t-t') = 0. (147) 

The corresponding diagrams are: 
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Eg(k,t-0 






For clarity, let us split the vertex made with W into two parts: 



The cross on a line with impulsion q stands for W(q) and the black triangle 
stands for k ■ q, k being the momentum of the outgoing arrow. Then we have 







and 



pw(k) 







Putting this altogether leads to 

£ Wp (k, t-t') = -(3 W(k) E M (k, t - t') 

for p = 2. 



(148) 
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Now let us assume that n > 4 is even and (|148|) holds for any even p less 
than n — 2. We have in frequency: 



n-2 



A (n) (k, u) = °o s(n " 9) ) (k, w) + (S w G A (n_9) ) (k, w), (149) 

g=2 

where we have set £( r ) = A( r ) = for odd r due to the absence of tadpoles. 
Diagrammatically: 



A )= E } E + r E A 








Let us consider U = A Q B, where A = £ (<?) and B = or £ = A( n_ «). 

We will show that U verifies (I148|) . which will be enough to show that A^ 
also verifies it due to (|149jl . We have: 

Up P (k,u) = Ap P (k,u) Ro(k,u) Bp P (k,u) 

Up P (k, to) = Atf(k, uj) R* (k, u) B pp (k, u) + A pp (k, u) C (k, u) B pP (k, to) 
+Ap P (k, lu) R* (k, u) Bp P (k, u). (150) 

However using the hypothesis made at order n — 2, A and B verify (|148|) . In 
addition, Q verifies (|141|) (note the absence of the minus sign compared to 
(|104|) ). Thus using (|141|) . one can express the terms with indices pp in the 
right hand side of ()150|) in term of those with indices pp and pp. Eliminating 
C in profit of G and Gq, fj!5(Jj) gives 



U (k w) - U PP (Kw) + Up P (Kw) 



Let us rewrite the above induction in such a way that the link with similar 
diagrammatic proofs in the case of additive noise is clearer (see [33] fc> r a 
diagrammatic proof of FDT in that case). The formal solution of (J149)) is: 

n— 2 r 

A?>,o;) = J2 E E £ } ( k ^) II (^(k,^) S^(k,o;)). (152) 

r=l iiH \-i r =n s=2 

Now we use the useful identity 

& (ZXZ2 {■■■)z n )= U{ Zl )4{- ■ ■ )z* n -Z 1 ^{z 2 ) Z*(- ■ ■ )z* n + Zl Z 2 £(* 3 ) (• • • )<- • ■ 

to get the following diagrammatic identity: 
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Again using the hypothesis made at order p < n — 2 and (jl41|) the result is 
obtained at order n. The importance of the minus sign in (J1U4)) compared to 
(ll(J3|) is clear here. As said above, this establishes (|1U4|) for class A^ n \ Thus, 
by subtracting diagrams of classes A^ n ' and from diagrams of Q^ n \ (jl04|) 
is established at order n. 




A. 2. Proof of FDT in the quadratic expansion of the 
free-energy 

The proof will follow the lines of the previous one. We shall proceed by 
induction. Here we notice that the diagrams of x contributing to FDT at 
order n are of order n — 1. So the first diagram corresponding to \ is for 
n = 2: 

x (2) = - 

Using the manipulations of the diagrams at order 2 made in the previous 
section, it is straightforward to prove FDT at this order. Then let us assume 
that FDT has been proved at any order p < n — 2. The SD reads for the 
diagrams of order n exactly: 

(iu + p Tk 2 W{k)) C (n) (k,cj) = 2Tk 2 po Q [n) {k,uj) 

+ E Hr^) 0«(k,a;) + 47 p) (k,^) C^(k,u)_ . 
p 

Thus, using the assumption made at order p < n — 2 to transform the self- 
energy part and (jl(J3|) to transform W(k) (k,u), we get 

- cW(k ^ =^,p(4r"(M) c«(k, U )) (is3) 

+ 2iT Po k 2 %g( n \k,Lj). 

Now we remark that we have 
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which shows that the term of Y* pp Q in (J153|) is nothing but the anomalous 
response %■ Thus the FDT is proved at order n. 



Appendix B. Causality with the extra fields 

In this appendix, we discuss briefly the difficulties in verifying that causality 
is verified by the SD equations in the presence of the fields 6 and 9. Causal- 
ity is mainly insured by the form of the bare propagator and the identities 
(|73|75f 75|76|) . However some of the diagram vanish due to the Ito prescrip- 
tion in a rather subtle way. We explicit this at the order of one-loop, the 
generalization to higher orders being straightforward. It is not difficult to 
be convinced that at one-loop the only diagram which may eventually cause 
some difficulties comes from the convolution (£ pp ■ C pp ) (k, r). This diagram, 
shown in Fig. |2] is proportional to 

J ^ k • (k - q) J™ dt dt [6(t)C>(q, t)} 6(-t)C>(q - k, -t)C pp (k, r-t). 

(154) 

However due to the Heaviside functions the integrals can be restricted to 
[— e, e], with e > 0. In addition 

J -j^L k - (k - q) ^ dt dt [6(t)C pe (q, t)] Q(-t)C pe (q - k, -t)C pp (k, r-t) 
« J ^k-(k-q) J" dtd t [Q(t)C pe (q,t)}C p e(q-k,0 + )C pp (Kr) (155) 

= J k • (k - q)C>(q - k, + ) C pp (k, r) [Q(t)C p9 (q, t)}^ 

= 0. 
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Figure 2: One- loop diagram contributing to S pp , which vanishing is the less 
obvious. 



Appendix C. Derivation of dynamical and static 
equations 

In this appendix we sketch the derivation of equations (|87|) and (j93j) as an 
example. Other dynamical or static equations can be obtain following the 
same routes. We start from the SD equation 

(G 1 -G-S-G) /3/3 (k,r)=5(r) (156) 

for any value of r. We have 

(E • G%p(k,T) = (E pp ■ C pp )(k,r) + (5> ■ C^)(k,r). (157) 

Indeed Cpp and C^- vanish by causality. We then get: 

(S ■ G)#(k, r) =^ | T <ft [<9 r E^(k, r - t)C>(k, f) 
+E /w (k,r-t)C7 M (k,t)]. 
Integrating by parts, one gets 



(E.%(k,r) = ^ |7) ^(k,i = 0) 
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T 

In addition we have 



/ dt [E^(k,T-O^C pe (k,t) + E^(k,T-f)C ee (k,t)]. 
Jo 



(Go 1 ■ G)^(k, r)5 T C^(k, r) + Po k 2 C ef) (k, r) 



Equating the terms proportional to 5(r) in (J156j) one gets 



158) 



(159) 



1 onk 2 (160) 

-9 r (0(r + )C p ,(k,r)) + ^C, e (k,r). 



Ic>(k,0) = l, (161) 
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and taking the limit r — > + , gives (|93|) . Finally (J87|) is obtained by taking 
t > 0. All other equations for the dynamical evolutions and the statics can be 
derived in the same way. When causality is not enough to restrict explicitly 
time integrals between and r, one can verify that in all cases FDT makes 
it possible to combine together different contributions of the same equations 
to finally end up with integrals between and r. 

Let us add here that careful analysis of the self-energies shows that E^ 
has a tadpole contribution. However this tadpole can be eliminated by adding 
a linear term —A J <i 3 x5p(x) to the entropic part of the free energy and A to 
the potential, with a suitable value of the constant A. 



Appendix D. Proof of the linear dependence 
of the Schwinger-Dyson equations 

As we have already mentioned, it may appear unnatural to describe the 
evolution of 3 correlators with 4 dynamical equations. A series expansions 
at low r > of these equations makes this clearer. When expanded in 
series, they become a cascade of equations for the successive derivatives of 
the correlators at zero time difference. We start by expanding at first order, 
from which it is easy to guess what is going on at higher orders. At order 1 
in t. (jHRjl . (j%7 jl . IftB jl and read respectively: 

C pp (k, 0) + p k 2 C>(k, 0) + r [c pp (k, 0) + p k 2 C>(k, 0) = 

r [£fc(k, + )C pp (k, 0) + Z p9 (k, + )C>(k, 0)] (162) 
C>(k, 0) + p k 2 C e ,(k, 0) + r [c>(k, 0) + p k 2 C^(k, 0) 
-£fc(k, 0+)C>(k, 0) + r [S^(k, 0)C ee (k, 0) - £ §e {k, + )C p9 (k, 0)} (163) 
W(k)C pp (k, 0) - C p9 (k, 0) + r [w(k)C pp (k, 0) - C p9 (k, 0) 

i£^(k, 0)C pp (k, 0) + r £^(k, 0+)C>(k, 0) (164) 



T 



W(k)C p9 (k, 0) - C ee (k, 0) + r W(k)C p9 (k, 0) - C flfl (k, 0) 



E fl - fl (k, + )C e ,(k, 0) - -E w (k, 0)C pe (k, 0) 



. (165) 
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In addition, the SD equations have an apparent singularity at r = which 
comes from the 8(r) in the RHS of ()54|). This gives an initial condition: 
C p g(k, + ) = T. Thus there are 5 equations at order 0, which fix the values 
of C pp (k,0), C pe (k,0), C ee (k,0), C pp (k,Q) and C p8 (k,0). At order 1, there 
are 4 equations but 3 quantities only to be determined, namely Cgg(k, 0), 
C pp (k, 0) and C p g(k, 0). Remark that the self-energies and their first deriva- 
tives appear in the equations. However, as they can be expressed in terms 
of the correlators, it can be checked that the successive derivatives of the 
self-energies can be expressed in terms of the quantities already computed 
at previous orders. This guarantees that at every order self-energies do not 
give extra variables to be determined. Now we show that in fact one of the 
equations obtained by identifying the terms of order r is trivially satisfied 
by the solution of the equations at order 0. We focus on the term propor- 
tional to r in the LHS of (jl64J) . We then express this term by using a linear 
combination of the terms of order of ()162|) and (jl63|) : 

W(k)C pp (k, 0) - C pe (k, 0) = - Po k 2 [W(k)C pp (k, 0) - C>(k, 0)] (166) 

+£^(k,0 + )C^(k,0). 

/.From order of ()164|) . the terms in brackets vanishes, and then we get: 

W(k)C pp (k,t) -C pe (k,t)£^(k,0 + )C>(k,0), (167) 

which corresponds to the terms proportional to r in (|164|) . Therefore the 
number of equations obtained at order r is equal to the number of variables 
to be determined at this order. 

The non perturbative generalization of the previous approach comes from 
the following remark: the SD equations form a linear system of equations 
which unknown variables are the correlators and coefficients are the compo- 
nents of Gq 1 and E. The solution of this system of equations is found easily 
using the Laplace transform. The SD equations read in Laplace transform: 



(168) 




^ (k ' Z) + (/?ok2 + ^P^)Coe(K z) = T (169) 
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±£ §§ (k,z)c pp (k,o+) 



T 



C p e(k, z) 



;i70) 



+ W{k) + 



«S^(k, z) - £^(k,0+) 



C pp (k,z) 



W(k) C p9 (k,z) 




T 



C p9 (k, z). 



For better clarity, we write formally this system as follows: 

zAC pp (k, z) + BC pe (k,z) 
zAC pe (k,z) 



- w ( k >°) 
T 



DC pp (k,0) 




EC pp (k, z) 
EC p6 (k, z) 



hBCg 6 (k, z) 
AC p9 (k,z) 
AC e9 (k,z). 



(171) 



(172) 
(173) 
(174) 
(175) 



The identity E (RHS) 1 - A (RHS) 2 - zA (RHS) 3 - B (RHS) 4 = 0, where 
(RHS)i stands for the RHS of the ith equation above, is trivially verified. It 
remains to prove that the LHS are linked by the same relation. Gathering 
the terms of E (LHS)! - A (LHS) 2 - zA (LHS) 3 - B (LHS) 4 = (with 
obvious notation), one gets: 



W(k)--E w (k,0 



C pp (k,0)=T. 



[176) 



This is precisely the static equation (JiMj). and the proof is complete. 



Appendix E. General fields transformation and 
field theory in the case of Brownian dynamics 
for the density field 

We have seen above (|4.1|) that one can make each of the two symmetries of the 
system linear independently. One can also render them linear simultaneously 
using the same method: introducing additional fields. 
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The starting point is the action (|64j) : the field transformation is given by 
Ti. We introduce 4 additional fields to linearize the symmetry U: 6 = SJ-'/Sp 
and rj — V • (pV^) as well as the conjugated fields 9 and r). The final 
transformation is: 



Px 
Ipx 
$x 

fx 

fx 

pj- 

0, 

fa 
fa 



Px + fx 

i>% + \dtPx 
^x + Tf x 

-f x + Tf x + ip x + Tp x 

fx 

—px + \®x 
~1px + ^Vx 

-4>x + Ox 

fjx + \®x -px~ \^x 

6x — \dtpx 

fx 
fx 



(177) 



:i78) 



The action equals to the integral over x, t of: 

-pxdtPx + i> x Q x ~ Tip x p x + f x (V ■ {p x Vf x ) + ^d t p x ) - p x V(j) x ■ V(j) 

M. 

'V' 



+M*Px - v ■ { Px Vp x )) + e x (e x -jf)+f] x ( Vx - v • (p,w,)) . (179) 



The structure of the two-fields correlators matrix for r > is schematically: 



T 




TC ff 

-Cff 


TCff 


°pf 

T 





n ^pf c pp 
^pp — 




Cpf 

Cpf 


Cpf 


Cpp 




/'/' 
T 




c pf 

T 

'Cpf 


°pf 

T 



Cpp 
T 





c pf 





TCff 

Cff 

c ff 

TCff 










Cpf 












C 



ff 











Cpf 



T 2 C ff 

-TCff 


T 2 Cff 


Cpf 




_C_pp_ 
T 




c pf 
. T 
c pf 
T 





C, 



p.! 



JUL 

T 



Cpp 















Cpf 





~Cpf 






















Cpp 


Crjfj 


Cryq 
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where the ordering of the fields is (p, p,ip,ip, f, f ,9,8,rj,rj). The fermionic 
fields and <p, which correlators with other fields vanish, have not been in- 
cluded. The role of these fermionic fields will be explained later on a practical 
example. For completeness, the bare propagator (still without the fermionic 
fields) is given in Fig. El Some of the correlators defined above are singular 



2k 2 p T 






iT 


2iw 


(1 







2k 4 p 2 TW 





2k 4 p 2 TW 


v; 2 +k 4 p' 2 W 2 


w-ik'pcW 


iw+k 2 p a W 


w+ikipcW 


w 2 +k 4 p 2 W 2 


iw+k 2 p a W 


w 2 +k 4 plW 2 


w 2 +k 4 plW 2 










W 













k 2 p W 





k 2 p W 


w+ik 2 p W 


-iw+k 2 p„W 


k 2 p T{w+ik 2 p W) 


-iw+k 2 p a W 


-iw+k 2 p a W 























k 4 p 2 W 





k 4 p 2 W 








w+ik 2 p W 


T(v,+ik 2 p W) 






-iw+k 2 p W 




-iw+k 2 p W 


iT 


w 

















ik 2 p TW 





k 2 p TW 


ui-ik 2 p W 


hu+k 2 p c W 


w-ik 2 p W 


k 2 p a (iw+k 2 p a W) 






w-ik 2 p W 


w-ik 2 p a W 




iw+k 2 p W 


2iw 








2w 2 


1 







2ik 2 p wW 


(] 


2ik 2 p wW 


w 2 +k°p'pv 2 


k 2 p T{iw+k*p W) 


T(w-ik* pvW) 


k 2 p (-ivi+k 2 pcW) 


k 2 pvT(w 2 +k 4 p 2 I W 2 ) 






T{w-ik*p W) 


w 2 +k 4 p 2 W 2 


in 2 +k 4 p 2 W 2 
















(P 











(] 


























-1 


























-1 










ik*p 2 W 


-iw+k 2 p a W 


-iw+k 2 p a W 


T(-iw+k 2 pvW) 


-iw+k 2 p a W 


w+ik 2 p W 


2k 4 p 2 TW 


k 2 pcW 


k'p 2 „W 


k 2 p TW 


2ik 2 p a wW 










2k e p%TW 2 


-1 


2k 6 p°TW 2 


w 2 +k*p 2 a W 2 


iw+k 2 p a W 


iw+k 2 pcW 


-iw+k 2 p a W 


w 2 +k 4 p' 2 W 2 


iw+k 2 p a W 


w 2 +k 4 p 2 W 2 


w 2 +k 4 p 2 W 2 


























-1 








2k*p 2 TW 


k 2 poW 


k A p\W 


k 2 p TW 


2ik 2 p a wW 








k 4 plW 


2k G plTW 2 





2k e p 2 a TW 2 


w 2 +k«p 2 W 2 


iui+k 2 p W 




-im+k 2 K W 


w 2 +k*p 2 W 2 






iu!+k 2 p W 


w 2 +k 4 p 2 W 2 




w 2 +k 4 p%W 2 




















(1 





-1 





-1 



Figure 3: Bare propagator obtained from the action f)179|) . where W stands 
for Wik). 



at r = 0. Indeed for r G M, we have C//(k, r) = a(k) 5(r) + (Tj^k, r) 
and Cjj = b(k) S(r), where C r ^ 9 is regular (i.e. continuous by parts) and 
even in time. We also have C^(k, r) = ^d T (6(r)C p /(k, r)) = — CWk, r), 
C^(k, r) = -C„{k, r) = c(k) 6(t) + TC r f f(k, r) and C vfj (k, r) = -5(r). 

We conclude this appendix with a discussion about causality. One can 
show using the way S transforms under {Ti,Ui} that S pp vanishes, which 
guarantees causality to hold non-perturbatively. Furthermore, it also hold 
perturbatively, as a consequence of the form of the bare propagator and of 
the Ito discretization. We first show how it holds at one-loop order, the 
generalization to higher orders being straightforward. 

At one-loop, S pp has three non-trivially vanishing contributions, shown in 
Fig. |3J The diagram involving the ghost loop is identical to the one involving 
the loop made with propagators Cjt, with opposite sign and with instead 

of C°.t. However the bare propagators C°- and C°~ are identical and the 

// <t>4> ff 

equations giving the renormalization of and C ^ are also identical, thus 
the contributions of both loops in the self-energy - which are divergent - 
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Figure 4: Diagrams contributing to S pp at one-loop. The fields involved in 
the Wick theorem are indicated around the vertices. 

cancel each other exactly. As usual in the Fadeev-Popov method, the role of 
the fermionic fields is to remove the volume of the "gauge" ensemble which 
is infinite. One can easily check that loops of ghost propagators annihilate 
at all orders with corresponding loops of propagators Cjf. Following the 
steps of appendix B, the remaining diagram in Fig. 0]may be shown to give 
a vanishing contribution when involved in the dynamic equations. 

Appendix F. Dynamical equations for Fluctu- 
ating Nonlinear Hydrodynamics 

In this appendix we give the derivation of the dynamic equations for fluc- 
tuating nonlinear hydrodynamics. The calculus and the ideas behind are 
the same as the corresponding for BDD although somewhat more cumber- 
some due to a larger number of fields. We start with the Schwinger-Dyson 
equations and use time-reversal to simplify them. 

This time (J66|) applied to the transformation U\ gives the following equa- 
tions for correlators: 



Cpp(k,r) 


T 


-C>(k,T) 


(181) 


CW(k,r) 


e(r) 

T 


'-C pv (k,r) 


(182) 


C p §(Kr) 


G(t) 
T 


-d T C pp (k,T) 


(183) 


C pil {k,T) 


Q(r] 
T 


-d T C pg (k,T) 


(184) 



(33 



Cgp(k, T 
C 99^ T 

C gi ,(k,r 

C 6p (k, r 
C 9g (k,T 

C d v(k,T 

C V p(k,r 
C vg (k,r 

C v {,(k, t 



T 
9(r 



T 
6(r 



T 
6(r 



6(r 



T 
0(r 



T 
6(r 



T 
6(r 



T 



0( 



r 



T 
6(r 



T 
6(r 



T 
0(r 



C ge (k,r) 
C gv (k,r) 
d T C gp {k, t) 
d T C gg (k,r) 

C ee (k,r) 
C ev (k,r) 
d r C 6p (k, t) 
d T C dg (k,T) 

C ve (k,r) 
C vv (k,r) 
d T C vp {k,r) 
d T C vg (k,r) 



and (fKUJ) yields the following identities for self-energies: 



s /5p(k, t) 
s /3 5 (k,r) 
S p e(k,r) 
S/3«(k, t) 



e(T)E rf (k,T) 



-a r [e(r)E / »(k,T)] 



T 

1 

T 



0(r)S pp (k,r) 
e(r)E pg (k,r) 



(185) 
(186) 
(187) 
(188) 

(189) 
(190) 
(191) 
(192) 

(193) 
(194) 
(195) 
(196) 

(197) 
(198) 
(199) 
(200) 
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s 3p( k > r ) 
S^(k,r) 



T 



T 



9 T [0(r)%,(k,r)] 
1 0(r)S^(k,r) 



T 

■ie(r)E M (k,T) 



(201) 
(202) 
(203) 
(204) 



S e >(k,r) 
^ 9 (k,r) 

s ee( k > r ) 
S fc (k,T) 



-e(r)«9 T £^(k,r) 

ie(r)a r E to (k,T) 

-ie(r)S^(k,r) 



0(r)E dd (k,r) 



(205) 
(206) 
(207) 
(208) 



E#p(k, t) 
Es g (k, t) 
S«e(k, r) 
Ej>i,(k, t) 

One can get some additional identities: 

S^(k,r) = 

^(k,r) = 

Efo(k,r) = 

Ej)p(k, r) = 

Efig(k, t) = 

^(k,r) = 



0(r)9 T S^(k,r) 



T 

i0(r)9 T S,. o (k,r) 
-i0(r)S^(k,r) 
-ie(r)^(k,r). 



s ra(k,r) 
-S^(k,r) 
-%(k,r) 
— E^(k, r) 
E^s(k, t) 
%(k,r), 



(209) 
(210) 
(211) 
(212) 



(213) 
(214) 
(215) 
(216) 
(217) 
(218) 
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and similar ones for correlators: 



C gp (k; t, s) 
C dp (k; t, s) 
Cg g (k;t,s) 
C vp (k; t, s) 
C vg (k;t, s) 
C ve (k;t, s) 



C P g(k;t,s) 
C p8 (k;t, s) 
C g0 (k;t,s) 
C pv (k;t,s) 
C gv (k; t, s) 
C fl „(k;t,s). 



(219) 
(220) 
(221) 
(222) 
(223) 
(224) 



All these identities reduce the number of independent correlators to ten, 
which are C P p, C P g, C P } C pv , Cgg ,Cg0 } C gv , Cgg, Cg v and C vv . 

In the case of FNH, there are in principle 64 Schwinger-Dyson equation. 
We write 16 of these equations, the other being trivially linear dependent on 
these: 



d T C pp (k, t) - ip kC vp (k, 



1 

T 



dtZ p§ (k,r-t)d t C pp (k,t) 



T 



- I dtEp(k, t - t)d t C gp (k, t)-- I dtEtf(k, t - t)C ep (k, t) 



- I dtX pg (k,r-t)C vp (k,t) 



(225) 



d T C pg {k, r) - ip Q kC vg (k, r) = - / dfE^(k, r - t)d t C pg (k, t) 

1 Jo 

+ i / dtEpik, r - t)d t C gg (k, t) - i / dtEppik, r - t)C 9g (k, t) 



TJ 



dtE pg (k, r - t)C vg (k,t) 



(226) 



d T Cp 9 (k, t) - ip kC v g(k, r) = E^k, r) 



S^(k, r - t)d t C p9 (k, t) + S pf ,(k, r - t)d t C g6 {k, t) 



+ T I dt 



£ w (k, r - t)C fle (k, t) - S p§ (k, r - t)C„ e (k, t) 



(227) 
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d T C pv (k, t) - ip kC vv (k, r) = E^(k, r) 

i / dt Z. § (k,T-t)d t C^(k,t) + Z^(k,T-t)d t C gv (k,t) 



- Epp(k, r - t)C 6v (k, t) - Ep-§(k, r - t)C vv (k, t) 



(228) 



d T C gp (k, r) - ip kC ep (k, r) + LC vp (k, t) + - dtZ § (k, r - t)d t C pp (k, t) 

1 Jo 

+ dtE g ,(k,r-t)d t C gp {k,t)-^ [ dtE gp (k,r-t)C 6p {k,t) (229) 



T 



<9 r C 99 (k, r) - z> kC eg (k, r) + LC vg (k, r) = +- / dtE-(k, r - t)<9 t C P5 (k, t) 



(230) 



T 

dtE^k, r - t)<9 4 C 99 (k, t) - 1 / dffi §(3 (k, r - *)C^(k, i) 



T 



TJ 



rftE § g(k,r - t)C„ ff (k,t) 



9 r C ge (k, r) - ip kC ee (k, r) + LC v0 (k, r) = E^(k, r) 
1 



dt 



E -(k,r - t)^C>(k,t) - E §fi (k,r - t)^C pfl (k,t) (231) 



£ w (k, r - t)C fle (k, t) + E^(k, r - t)C„ fl (k, t) 



(232) 



+ 



<9 T C 5U (k, r) - ip kCov(k, r) + LC w (k, r) = E §fi (k, r) 

i / dt E^(k,r-09 t C /w (k,0 + Ep(k,T-t)a t C fl „(k,f) (233) 



- Eg«(k, r - t)C e „(k, t) - S^(k, r - t)C w „(k, t) 
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C ep (k, r) - W(k)C pp (k, t) = -- p w (k, 0)C pp (k, r) + E^k, 0)C gp (k, r)] 



T 



dt 



E fll (k, r - t)<9 4 C pp (k, t) + E fe (k, r - t)<9 4 C 9P (k, t) 



(234) 



E^(k, r - t)C*„(k, t) - E^(k, r - t)C wp (k, t) 



C e9 (k, r) - W(k)C P9 (k, r) = -- p w (k, 0)C P9 (k, r) + E fe (k, 0)C g9 (k, r)] 



eft 



E w (k, r - t)9 t C pg (k, t) + E fe (k, r - t)ftC OT (k, t) 



E 0> (k, r - t)C eg (k, t) - E^(k, r - t)C^(k, t) 
C ee (k,r) - ^(k)C>(k,r) = -E^t) 



- [E w (k, 0)C>(k, r) + E fc (k, 0)C>(k, r)] 



(235) 



(236) 



rj dt 



E §§ (k, r - t)d t C pe (k, t) + E fe (k, r - t)0 t C flfl (k, t) 
E e> (k, r - t)C M (k, *) - E fe (k, r - *)0*(k, *) 



C to (k, r) - ^(k)C p ,(k, r)-E fe (k, r) 



(237) 



1 

T 
1 



- [E w (k, 0)C p ,(k, r) + Efc(k, 0)C fl ,(k, r)] 

E w (k, r - t)^C p „(k, t) + E fe (k, r - t^C^k, t) 

E^k, r - OOUk, t) - E^(k, r - t)C vv (k, t) 



C vp (k, r) - -C 9P (k, r) = ~ [E^(k, 0)C pp (k, r) + E M (k, 0)C 9P (k, r)\ 
Po J- 



T 



dt 



E-(k, r - t)^C pp (k, t) + E co (k, r - t)9 t C flp (k, t) 



(238) 



Ej) p (k, r - t)C ep {k, t) - E^(k, r - t)C vp (k, t) 
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C vg (k,r) - -C gg (k,r) = -I[S^(k,0)C w (k,r) + E cc (k, 0)C gg (k, r)] 
Po J 

+ ^ f dt E^(k, r - t)^C w (k, t) + E e{i (k, r - *)<9 t C gg (k, t) (239) 
J ^o L 

- E C(S (k, r - t)C 0p (k, t) - Eyg(k, r - t)C vg (k, t) 



+ 



C ve (k,r)--C ge (k,r) + ^ d§ (k,r) 
Po 1 

i [s^(k, 0)C^(k, r) + E M (k, 0)C^(k, r) 

1 



(240) 



E M (k, r - t)d t C pd (k, t) - E M (k, r - t)0 t C flfl (k, t) 



- E^(k, r - t)C flp (k, t) + E fi§ (k, r - t)C„ fl (k, t) 



C m (k, t) - — C gv (k, t) + — E M (k, r) 
Po T 

^ [E^(k, 0)C p ,(k, r) + E^(k, 0)C ff ,(k, r) 



(241) 



+ y / dt 



E^(k, r - t)^C p „(k, t) - E«H)(k, r - t)$C fl „(k, t) 



E C/3 (k, r - t)Cg p (k, t) + E, g(k, r - t)C vv (k, t) 



As for BDD, the number of independent correlators is smaller than the num- 
ber of equations, and there are here 6 redundant equations. The extension 
of the proof of appendix D to the present equations is straightforward but 
very painful. 

colorblue These are the exact non-perturbative dynamical equations pre- 
serving FDT. One can then use different approximation schemes for self- 
energies to concretese the equations. It is worth noting that what ever the 
approximation is the FDT is alwas verified due to theway the equations were 
derived. 



69 



References 

[I] C. A. Angell, Science 267, 1924 (1995). 

[2] Recent reviews: P.G. De Benedetti and F.H. Stillinger 
Nature 410, (2001) 267; 

M.A.Ediger Annu. Rev. Phys. Chem. 51 (2000) 99. 
[3] S. P. Das Rev. Mod. Phys. 76 785 (2004). 

[4] D. R. Reichman and P. Charbonneau, J. Stat. Mech. (2005) P05013. 

[5] W. Kob, in "Slow relaxations and nonequilibrium dynamics in condensed 
matter" , vol. Session LXXVII of Les Houches Summer School (ed. J.-L. 
Barrat, M. Feigelman and J. Kurchan), published by EDP Sciences and 
Springer. 

[6] U. Bengtzelius, W. Gotze and A. Sjolander, J. Phys. C 17, 5915 (1984). 

[7] E. Leutheusser, Phys. Rev. A 29, 2765 (1984). 

[8] W. Gotze and A. Sjolander, Z. Phys. B 65 415 (1987). 

[9] S. P. Das and G. F. Mazenko, Phys. Rev. A 34, 2265 (1986); S. P. Das 
Phys. Rev. A 42, 6116 (1990). 

[10] J.-P. Bouchaud, L.F. Cugliandolo, J. Kurchan and M. Mezard, in "Re- 
cent progress in random magnets", A. P. Young ed., World Scientific 
(1997). 

[II] J.-P. Bouchaud, L.F. Cugliandolo, J. Kurchan and M. Mezard, Physica 
A 226, 243 (1996). 

[12] L. Berthier, J.-P. Garrahan, Phys. Rev. E 68, 041201 (2003). 

[13] K. Miyazaki and D. R. Reichman, J. Phys. A: Math. Gen. 38, L343-L355 
(2005). 

[14] T.R. Kirkpatrick and D. Thirumalai Phys. Rev. A 37, 4439 (1988). 

[15] S. Franz, G. Parisi, J. Phys.:Condens. Matter 12, 6335 (2000). C. Donati, 
S. Franz, S.C. Glotzer, G. Parisi, J. Non-Cryst. Sol., 307, 215-224 (2002). 



70 



[16 
[17 

[18 
[19 

[20 
[21 

[22 
[23 
[24 

[25 
[26 

[27 
[28 
[29 
[30 

[31 
[32 



G. Biroli and J.-P. Bouchaud, Europhys. Lett. 67 (2004) 21. 

C. Toninelli, M. Wyart, L. Berthier, G. Biroli, J.-P. Bouchaud, Phys. 
Rev. E 71, 041505 (2005). 

H. C. Andersen, PNAS 102 6686 (2005). 

A. Latz, J. Phys: Condens. Matter 12, 6353 (2000); A. Latz, J. Stat. 
Phys. 109, 607 (2002). 

M. Fuchs ad M. E. Cates, Phys. Rev. Lett. 89, 248304 (2002). 

K. Miyazaki and D. R. Reichman, Phys. Rev. £66 050501 (R) (2002); K. 
Miyazaki, D. R. Reichman and R. Yamamoto, Phys. Rev. E 70, 011501 
(2004). 

C. Jarzynski, Phys. Rev. Lett. 78 2690 (1997). 
G.E. Crooks, J. Stat. Phys. 90 1481 (1998). 

P. M. Chaikin and T. C. Lubensky, "The Principles of Condensed Matter 
Physics", Cambridge University Press (2000). 

D. S. Dean, J. Phys. A 29, L613 (1996). 

P. N. Pusey, in "Liquids, Freezing and Glass Transition" , J. P. Hansen, 
D. Levesque, J. Zinn- Justin, Eds. (North-Holland, Amsterdam, 1991), 
pp. 763-942. 

K. Kawasaki and T. Koga, Physica A 201, (115). 

. C. Chamon and L.F. Cugliandolo, unpublished, 2004. 

T. V. Ramakrishnan and M. Youssouff, Phys. Rev. B 19, 2775 (1979). 

J. Zinn- Justin, "Quantum field theory and critical phenomena" (Claren- 
don Press, Oxford, 1989). 

P. C. Martin, E. D. Siggia and H. A. Rose, Phys. Rev. A 8, 423 (1973); 
H. K. Janssen, Z. Phys. B 23, 377 (1976); C. De Dominicis, J. Phys. 
(France), Colloq. 1, C-247 (1976). 

C. De Dominicis and P. Martin, J. Math. Phys. 5, 14 and 31 (1964). 



71 



[33] S. K. Ma, "Statistical mechanics", world scientific (Singapore, 1985).; 
S. K. Ma, "Modern Theory of Critical Phenomena", Perseus Publishing 
(2000). 

[34] K. Kawasaki, Ann. Phys. 61, 1 (1970). 

[35] M. E. Cates and S. Ramaswamy, lcond-mat/0511260| (2005). 

[36] R. Schmitz, J. W. Dufty, and P. De Phys. Rev. Lett. 71, 2066 (1993). 

[37] W. Gotze and L. Sjogren, J. Phys.: Condens. Matter 1, 4203 (1989). 

[38] W. Gotze and L. Sjogren, Rep. Prog. Phys. 55, 241 (1992). 

[39] K . Kawasaki and S. Miyazima, Z. Phys. B: Condens. Matter 103, 423 
(1997). 

[40] G. Szamel, and E. Flenner. Europhys. Lett. 67, 779 (2004). 
[41] H. Osada, Probab. Theory Relat. Fields 112 53 (1998). 
[42] G. Biroli and A. Crisanti, unpublished. 

[43] J. P. Hansen, I. R. McDonald, "Theory of Simple Liquids", second edi- 
tion, Academic Press, London (1986). 

[44] T. Morita and K. Hiroike, Prog. Theor. Phys. 25, 537 (1961). 

[45] T. R. Kirkpatrick, P. G. Wolynes, Phys. Rev. A 35, 3072 (1987). 

[46] J. Berges, N. Tetradis and C. Wetterich, Phys. Rep. 363, 223 (2002). 

[47] P. Mayer, K. Miyazaki, D. R. Reichman, Cooperativity Beyond Caging: 
Generalized Mode Coupling Theory, cond-mat/0602248, 



72 



